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FOREWORD 


I 


Project Whirlwind 


Project Whirlwind at the Massachusetts Institute of Technology Diffita 
" Th ° r M 0r \ i ‘ * p0 , n, ° r ' d b * ,ht <*««« »f Naval Research undcrCon.rac 
diB^lal 4 co Jnu.a 0 j f C ’ Project are (1) the application of an electronic 

ln 8 matheT? of >rge capacity and very high speed (Whirlwind 1) to problem- 
^dd^l ’ V l ' nc *' engineering, simulation, and control, and (2) the studi 

and development of component reliability in Whirlwind I. 1 ' * 


The Whirlwind 1 Computer 




‘ . me nign-speed electronic digital type, in which ouantiti 

are represented as discrete numbers, and complexproblem, arc !Solved bit 

W„‘a; y ( 2 LX g : eVtoVe^a 1 ^"'* ITV" -moryf’l^wh" 

magnetic core, ‘'° r " ia! ‘ 0nf of ,wo '•‘rations of magnetic flux within ferr. 


1. QUARTERLY REVIEW AND ABSTRACT 


Jau.ing the past quarter 4b problems made use of the computer time allocated to 
the Scientific and Engineering Computation (SfcEC) Group. A procedure was developed to 
eliminate most of the need for manual intervention in reading in the tapes for different 
problems during a computing period, new routines for automatic curve plotting were intro¬ 
duced into the comprehensive system of service routines. 

Two simulated computers were developed, a three-address computer (TAC) and a 
single-address computer (SAC), which is a modification of the Summer Session Computer 
described in Summary Report No. 35. 

Work on the central computer by the systems group included modernization of the 
arithmetic-element control, reorganization of the power-supply control for a large portion 
of the terminal equipment, and equipment changes to facilitate trouble shooting. The addi¬ 
tion of two drum groups of auxiliary storage has increased the capacity from about 31,000 
to about 36,000 registers. 

Life tests involving tungsten-nickel Ail Cathaloy cathode sleeves, being run on 
5687 s and some special GE Z-2177's, show that there is more interface impedance present 
on the conducting sections than on the cutoff sections. This is the reverse of experience 
on active alloys in which silicon is the impurity. 

Of the 35 special programs given by MIT during the summer of 1954, two were 


d.Bit.) W TW. W 1«1? 1 !b US " nU T b ' ri . 0 ' 16 binary digit, (equivalent to about 5 decimt 
0 ..JL 1 V A g wa5 * rl ' c, * d *0 limit the machine to a practical ,ixe hu? 

sr. ; =i 

speed of the computer is 40 000 sinal** ' nary digits. Preser 

to about 20.000 multiDlication. „"« lf - add "“ °P'™‘'°n, per ,econd. equivaler 


directly concerned with digital cornputeition. and several others were of interest in the 
field of computer application. These courses, of 1 or i weeks duration, were largely for 
the benefit of people from business and industry. 
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2. MATHEMATICS, CODING, AND APPLICATIONS 


2.1 Introduction 

During the periodcovered by . hi, report 46 problem, made u.eot the computer time 
otted to the Scientific ard Engineering Computation (SfcEC) Group. Progress report, as 
submitted by the various programmer, are presented in numerical order in Section 2 2 
Of these, eighteen (185. 1,0. 1,3. 1,4. 1,5. „ 6 . 1,7. „ 8 . i0 0. 201. 202. 204 20 5 ‘ 

iOb 208. 10,. and 211, represent new problems .ha, are being de.c ribed for the firs, time! 
Thirteen problems (147. 14,. ,6,. 174. 175. ,85. „0. „2. „6. „7. 1,8. 202, and 20,) 

have been completed. One problem (143). previously reported a, completed, actually 
continued into the present quarter. 

Two simulated computers were developed for use in the 2-week summer session 

:z:Z'r< Com r;; : App “ ca,,on *-" Th *—* ««•—«.. c 0mP :;:: 

called TAG. is described in Section 5. 2 of this report. The other, a single -address com 
ron5?of d S SAC ' *' ‘ m0d ‘ f ‘ Cat,0n ° f the Summ " (SS) Computer described Sec- 

r”■■*- - »—-■ 

A procedure has been developed to eliminate most of the need for manual interven¬ 
tion in reading in the tapes for different problems during a computing period The proce 
dure ..described under Problem 100. 1, makes use of a specially prepared "director" tape 

."r-" 

Problem 100. ‘ rouUnes also be found under 


2.2 Problems Being Solved 

100 COMPREHENSIVE SYSTEM OF SERVICE ROUTINES 

- r ;••• * - -.* 

.cm used in CS U. the follow^ iTiZZ^on ■T ZlY'lZll "" 

ence. e lor thc feeder's conveni- 
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(m.n) number* shall mean number* which are of the form * * x • Z* where x i*an 
m-bina ry-digit number and y is an n-binary-digit number. For example, (24.6) signifies 
a two-register floating-point system dealing with numbers of 24 significant binary digits 
(roughly seven decimal digits) with magnitudes between 2^* and Z“^. 

Arithmetic involving these (m.n) numbers is carried out by means of (m.n) inter- 
pretive subroutines These subroutines enable the programmer to write coded programs 
using (m.n) numbers as easily as. or even more easily than, he might write programs in 
the single-length fixed-point (15,0) number system which is built into Whirlwind I. 

Director Tapes 

During normal operation of the comprehensive system of utility programs a pro¬ 
grammer submits a performance request for each groupof programs to be run on the com¬ 
puter. The performance request contains a complete description of the run. 

A set of standard instructions for use on the performance request has beendefined. 
Included among these are the following: 

e, Press the erase button 

fbl00-0-0 ( Place the designated tape in the input device (generally the photo¬ 
electric tape reader). 

ri. Press the readin button. (This reads in a tape placed in the input 

device and erases core memo.y if the erase button has previously 
been pushed.) 

The tapes involved in a performance request are prepared and assembled by the 
Tape Preparation Room. Actual computer operation is performed by trained operators in 
accordance with the instructions contained on the performance request. 

A partial mechanization of this process has been achieved with the introduction of 
director tapes. A director tape for a particular computer run is a punched paper Jape 
which is obtained by typing the performance request for the run on a Flexowriter. 

The comprehensive system of utility routines has been modified so that the major¬ 
ity of computer runs can be made under the control of a director tape. Initially the direc¬ 
tor tape and the run tapes (spliced together in the proper order) will be placed in separate 
input devices. The complete run will require only a single pushing of the readin button. 

The introduction of director tapes should lead to a more efficient use of computer 
time with less possibility of operator error. 

Additional Automatic Output Features 

Ten routines (five interpreted and five noninterpreted) have been written for use in 
automatic curve plotting. 

Two different "formats" are available for prescribing the ranges of the variables to 
be plotted: (i)FOC 1, for use when the range of variables involved is symmetric, and 
(i)FOC 2, for use when the range of variables is listed as program parameters. 
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Table 2-II. Current Problems Arranged According to the Mathematics Involved 
(* MIT Project on Machine Methods of Computation) 
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Point-by -point curve plc’ting is done by using ihe instruction (l)SOC si. The ab¬ 
scissa of the point is contained in register »l. and the ordinate of the point is assumed to 
be contained in the AC (or the MRA). 

Calibrated axes can be plotted on the scope by using the instruction (i)SCX. The co¬ 
ordinate. of Ihe origin and the calibration increments are listed a. program parameters. 
Uncalibrated axe, can be plotted on the scope by using the instruction (i)SUX. The coordi- 
nate» of Ihe origin are listed as program parameters. 

106 MIT SEISMIC PROJECT 

As discussed in various previous reports. Problem 106 is concerned with the in¬ 
vestigation of the use of statistical-analysis techniques in seismic-record interpretations 
and in particular in the separation of ■■reflection." from background interference on these 
records. A more complete desc r.pt,on of the problem and the approaches used is contained 
in Detection of Reflections on Seismic Records by Linear Operators" (Wadsworth. Robin¬ 
son, Bryan, and Hurley.-Geophysics. Vol. 18, No. 3 , July 1953 ). 

Because of the multiplicity of directions in which research on the above problem ha, 
gon- recently, we shall describe the computational aspects of Problem 106 in two ways 
first as they relate to these directions and secondly according to the mathematical tech¬ 
nique, they represent. We shall then describe recent results and finally indicate the prob- 
able course of future research. 

According to the type, of research they represent, our computation, fall into the 
following categories; 

1. Computation, designed to test criteria for automatic selection of optimum linear 
operators or linear filters; 

*’ ^' ompu,a «* on * designed totes, technique, for probing the statistical and frequen- 
cy behavior of a given input seismogram; 

3. Computations aiding the development of the relationships between discrete lin¬ 
ear operators and linear electric filters; 

9. Computation, designed to simulate and study various statistical noise model. 

5. Computation, for relating result, from 1 and 2 with geologic structure 

The type, of computation, involved in these categories are varied and overlap con¬ 
siderably. They include; formation of auto- and c ro, .correlation ma, rice ,. ma.rixmulti- 
plication, e.genroo. and eigenvector determination, discrete convolution, formation of auto- 
and crosscorrelations by several different formulas, power-spectra, estimates and ampli¬ 
tude and phase estimate, according to various formulas, a variety of linear and quadratic 
mixing techniques, numeric., differentiation and integration, frequency-distribution com- 
putation, series generation by moving summation, and others. 
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Some of the specific results accomplished during the past quarterly period includc 
the successful testing and use of new criteria mentioned under I, above, the development 
and successful testing of new probing ter hniqur s unde r above; the initiation of a noise- 
generation and study project, and the beginning of di ret t *.t rue tural intc rprctation programs. 

Our plans for the future are to pursue research a. indicated under all five catego¬ 
ries above, probably putting special emphasis on the linear operator - linear filter rela¬ 
tionship and statistical noise models, categories 3 and 4. 

107 (a) AUTOCORRELATION AND (b) FOURIER TRANSFORM 

INTEGRAL EVALUATION 

Routines developed under this problem by D. T. Ross of the MIT Servomechanisms 
Laboratory have been used by F. Raichlen of the MIT Hydrodynamics Laboratory in'the in¬ 
vestigation of tu rbulem-veloci ty flue’nation* in open-channel flow mrasurrH by means of a 
Pitot-tube pressure-cell combination. 

Previously the autocorrelation curves obtained from the Digital Computer Labora¬ 
tory had a delay-time spacing of 0.0015 second. This spacing didn't give sufficient infor¬ 
mation in the neighborhood of T =0. Therefore 1000 points were taken from these iccords 
with a time spacing of 0.00030J second, thus giving more information about the higher 
frequencies present in the turbulent-velocity fluctuations. 

The autocorrelation curves and mean-intensity spectra previously obtained were an¬ 
alyzed to determine the longitudinal scale of turbulence, L^. a length which is associated 
with the average eddy size. 

Since L^ is defined as L^ E U J* R(r)dT, 

where U = mean local velocity, and 
R(r) = correlation function, 

it would seem that the easiest way to ar rive at these values would simply be to lake the area 
under the autocorrelation curve between the curve and the line of zero correlation. Be¬ 
cause of the natural frequency of the Pitot-tube pressure-cell combination, the autocorrela¬ 
tion curve received from the Digital Computer Laboratory had periodicities superimposed 
on the autocorrelation curves for the random fluctuations. It would have been very difficult 
to average out these periodicities, with good precision, to arrive at the random part of the 
curve and hence the area between the random curve and the line of zero correlation. There¬ 
fore, an alternate method was devised whereby the theoretical curve for the mean-inten¬ 
sity spectrum, 

4L /U 

w(0 =---j-. 

1 \ (2W(L /Up 


II 
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was obtained from assuming the form of the autocorrelation curve, R(T) = .- T U/L x By 
fitting the former curve to the experimentally delermmed spectra, a series of values of 
L r and L^ (since UL r ) for different depths was determined. 

The problem then was to determine values of L x directly from the autocorrelation 
curves to check the L x values obtained from the mean-intensity spectra. Basically this 
was done by realising that the area under the mean-intensity spectrum and the r = 0 ordi- 
nate on the autocorrelation curve were the mean square of the turbulent-velocity fluctua¬ 
tions. By dele rmining the percentage of the area in the power spectrum representing pe¬ 
riodicities and reducing the r = 0 ordinate of the autocorrelation curve by this percent¬ 
age, a starting point for the theoretical autocorrelation curve was obtained. Values of L 
and hence L x were determined by fitting the curve e- TU/X x to the data, using as the z-rl 
or-hna.e the one explained previously. The results obtained have been found to agree very 
well with results of other investigators. 

120 T EM^PE R A TURE^V^ELCR~I Vy^^\S^ST^R*e*aM ^ INJECT ‘°» ‘ NT ° »‘ GH - 

This problem ,s connected with the development of a potential gas-turbine compo¬ 
nent. called an "ae rothe rmopressor, " which a net rise in stagnation pressure of a hot 
gas stream ,s brought about by the evaporation of liquid (water) injected intoa h.gh-velocity 
region of the flow. The concepts underlying the operation of the aerothermopressor are an 
outgrow,hof comparatively recent work in the field of gas dynamics, and it, intended func¬ 
tion the gas-turbine cycle is analogous to that of the condenser in a steam power plant. 

urther descriptions of this device may be found in earlier reports, beginning with Sum- 
mary Report No. 32, Fourth Quarter 1952. 

Because Of the complexity of the thermodynamic and dynamic processes whichoccur 
witnin the aerothermopressor involving the simultaneous acceleration and evaporation of 
liquid carried by an accelerating compressible How, any meaningful theoretical analysis of 
the problem currently being treated by Whirlwind I involves the simultaneous solution of 
seven nonlinear, ordinal y differential equations. 

During the past 15 months, a simple Euler method type of numerical integration was 
used satisfactorily forcompu.ing the performance of constant-a rea aero.he rmopressors un¬ 
der various operating conditions. Recently, however, U was found that such simple proce¬ 
dures were entirely inadequate for treating aero,hermopre.sor.of varying cro..-..c.io„,l 

area, and the resulting truncation errors led to results which were difficult to interpret 
The work carried out during the pas, quarter was devoted almost entirely to ihe develop¬ 
ment of a more eUborate numerical procedure which would eliminate serious trunca.il 
error and at the same time allow the useo, reasonable increment sizes so that the compu- 
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tation time would not be prohibitive. 

Two basic methods were studied (1) a forward-integration scheme involving back¬ 
ward differences in the fi-st derivatives of the dependent variables, and (2) the fourth- 
order Kunge-Kutta method. A third scheme, involving forward and successive integra¬ 
tion, was held in abeyant e pending results from the aforementioned method ,. Detailed 
studies were made with various increment sir.es for a typical case, and it was found that 
the more elaborate procedures did indeed eliminate the aberrations brought about by trun¬ 
cation error with the Euler method. Ihe Runge-Kutta method proved to be superior from 
points of view of speed, accuracy, and convenience, since no special starting procedure is 
required. For an accuracy of about 0. 3% in the pressure, the Runge-Kutta method was 
about six times as fast as the Euler method. 

The importance of the vanable-area aerothermopressor process has made it man¬ 
datory to incorporate the more accurate numerical procedure into the analysis. A new pro¬ 
gram is accordingly being written which, in addition to using the Runge-Kutta method, will 
exploit all the latest facilities of tlie comprehensive conversion system of Whirlwind l. Since 
the writing of this program will take considerable time, it is planned to continue computa¬ 
tions of aerothermopressor performance to a limited extent using the composite program 
which was used for testing the various numerical procedures. This program suffers pri¬ 
marily from inefficiency, since its basic logic was never intended for this purpose, nor was 
the auxiliary magnetic drum available at the lime it was written. 

The ae rothe rmopressor development program is beingcarried out at MIT under the 
sponsor ship of the Office of Naval Research and is directed by Professor Ascher H. Shapiro 
of the Department of Mechanical Engineering. The theoretical aspects of the problem treated 
by Whirlwind I are being carried out by Dr. Bruce D. Gavril. 

122 COULOMB WAVE FUNCTIONS 

Work was continued by A. Temkin of the MIT Physics Department on the Coulomb 
wave functions making use of their integral representation. It will be recalled that the reg¬ 
ular solution is given by: 
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A six-point recursion formula is being used (or the integration, and our solution 
was tested against some known results. We found that the accuracy increased as we de¬ 
creased the mesh size, but that extending the rang” of integration beyond a reasonable point 
did practically nothing. Some other con• .derations alsogive additional accuracy. Namely, 
near the origin the function tanh £ varies roughly as £ , whereas beyond £ * i, tanh £ is 
pract'cally constant. Thus the frequency of oscillations of the integrand around the origin 
is determined by | p -It) |. whereas further out it is determined only by 1 7 ). Since the 
mesh size should be inversely proportional to the frequency, we test to see if \l7)-p\ > irj, 
and if it is we use the mesh size corresponding to it in the region £ =0-^1. 

We are now testing for various values of the parameters to make sure that the rou¬ 
tine can handle them all. When this is done, we shall code the irregular solution and the 
first derivatives of both. We hope to be done with the coding by the end of the next quarter. 

123 EARTH RESISTIVITY INTERPRETATION 

This problem is concerned with a means of analyzing ea rih-resistivity data. The 
procedure is to fit a theoretically derived Slichter kernel to an observed kernel by system¬ 
atically relaxing the values of resistivities and thicknesses used in the theoretical kernel 
(see Summary Report No. 27). 

The program for this analysis is being modified in an effort to find, empirically, 
some optimum sequence of parameter-change multipliers. It was found that, because of an 
error in the program last used, the value, of these parameter, were increased only once, 
from 0.05 to 0. 1, after the 20lh iteration. In addition, it was found that the calculation of 
exponential, was not a, accurate a, desired for large arguments, and this is being im- 
proved. 

Several more kernel, have been run. The usual behavior has been convergence for 
a short while, and then rapid divergence, although one example converged rapidly and con- 
tinued making small oscillations about the proper solution. 

A hand check of the final result of the 91-step analysis described in Summary Re¬ 
port No. 38 indicate, that there i, some question as to the validity of the final result, since 
the actual error is greater than that calculated by WWI and is. in fact, slightly larger than 
the sup Which gave the correct solution. The reason for this lie, in the program and ha, 
not been positively located as yet. 

This work is being carried out by K. Vozoff of the Department of Geology and Geo¬ 
physics under the sponsorship of the Project for Machine Methods of Computation and Nu- 
merical Analysis. 



_ 
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126 A DATA-REDUCTION PROGRAM 

Problem 126 is avery large data-reduction programfor use in the Servomechanisms 
Laboratory. The over-all problem is composed of many component sections which will be 
developed separately and then combined at alaterdate. Thus far, efforts have been focused 
or the development of utility-type programs. These programs, which have been described 
in previousQuarterly Reports,include a fully automatic program to fit polynomials to arbi¬ 
trary empirical functions, a general-purpose Lagrange interpolation program, and a flex¬ 
ible and fairly elaborate post-mortem routine. The programs arc being developed by Doug¬ 
las T. Ross and William M. Wolf with the assistance of Miss Dorothy A. Hamilton, Servo¬ 
mechanisms Laboratory staff members. 

Three versions of the complete data-reduction program, each version an extension 
of the previous one, have been tested during this period. Most of the tests have been run 
using the mistake diagnosis routine (MDR) to extract over 30 intermediate results for ac¬ 
curacy and logical checks. Since the optimum settings for certain parameter. are unknown, 
the facilities of the MDR arc nearly indispensable, and the routine use of the MDR has aided 
the development of these programs considerably. The testing of the third version, called 
the basic data-reduction program, is nearly complete. When this basic version is com¬ 
plete the first two developmental versions will be discarded, and further elaborations of the 
basic version will be written. 

Work is progressing on another phase of this problem, using auxiliary in-out equip¬ 
ment both as an integral part of th** data-reduction programs for monitoring purposes and 
as an aid to experimenting with improved computational techniques. Intervention registers 
and scopes will constitute the working parts of this system with efforts being made to pro¬ 
vide the flexibility normally associated only with analog computers. Present work is con¬ 
cerned with the development of routines to decode intervention-register inputs and to util¬ 
ize various specialized scope outputs. After merging these routines with the data-reduction 
programs, an elaborate logging scheme, using magnetic-tape output, will be written for 
permanent records of all operations and interventions during an experimental run. 

131 SPECIAL PROBLEMS (STAFF TRAINING, DEMONSTRATIONS, ETC.) 

This problem has been set aside to account for the WWI time expended in demon¬ 
strating the computer to visitors at the Digital Computer Laboratory, in developing rou¬ 
tines to be used in these demonstrations, and in testing routines written as part of the train¬ 
ing of new staff members. 

During the past 3 months, 15 groups visited the Laboratory. The affiliation of 
some of the larger groups are given in Section 5. 3. 
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132 SUBROUTINES FOR THE NUMERICALLY CONTROLLED MILLING MACHINE 

During the first part of this quarter, debugging of a program for preparing milling- 
machine tape for airfoil templates was completed by J. H. Runyon of the MIT Servomech¬ 
anisms Laboratory. Much of the remainder of the programming and computer time was 
devoted to the development of programs for checking milling-machine tapes. These pro¬ 
grams are based on a library subroutine for reading, decoding and detecting errors in the 
tapes. The programs arc for obtaining printouts of decoded milling-machine orders, cum¬ 
ulative totals of orders, and indications of errors. Several tapes containing, in all, about 
7500 milling-machine blocks have been successfully checked. 


Ill 


SfcEC SUBROUTINE STUDY 


Several subroutines in the library have been replaced or modified. 

These are: 

LSR FU 1 Exponential Subroutine. New version is 55 registers long and uses one double- 
length or one buffer register for temporary storage. 1 The first eight terms of 
the power series for exp x are used lo evaluate the function. The subroutine 
can be used with any (3U-j,j) PA. 

LSR FU 2 Square-Root Subroutine. Two versions are now in the library. The first, LSR 
FU 2a, is 34 registers long and uses two double-length or two buffer registers 
for temporary storage. The square root is formed by Newton's method. It¬ 
eration continues until the relative error between the last two iterates is less 
than 2 . This program uses (24,6) arithmetic. The second subroutine, LSR 

$U 2b, is 27 registers long, with the same temporary storage requirement as 
FU 2a. In this subroutine, the square root is also found by Newton's method, 
going through four iterations for all positive numbers. This subroutine can be 
used with any (30-j, j) PA. In both subroutines, if x is +0, control is transfer red 
back to main program with x in MRA. If x is less than -0, the computer stops. 

LSR FU 5 Hyperbolic Functions. This subroutine calculates sinh x or cosh x according 
to where entered. The new version is 89 registers long and uses one double- 
length or one buffer register for temporary storage. 1 

LSR MA 4 Symmetric Matrix Diagonalization. This is a shorter and f*«t* 
program described in Summary Report No. 35, pp. 19-21. 


version of the 


Specified by preset parameter. 


The auxiliary buffer routine (SP 1) has been discarded, since it is nolonger useful. 
Two subroutines hav< been added to the library 


LSR SP 3 


and a divisor cFd- 
this i online calcul ates and stores the first two terms of the binomial serics 


30 x 10 Divide. With a dividend of the form a ♦ b • 2 


LSR SP 4 


(i) . i( b + i( - “ + !*)Z' i0 . 

'c' C 1 C ' C C c 2 

Because of roundoff in the subroutine and truncation in the series, there will be 
a maximum error affecting the last three binary digits, with a normal error of 
two digits being affected. 

30 x 30 Divide. The method of the subroutine is roughly that of the $S divide 
where the magnitude of the divisor is multiplied by increments, and the results 
are subtracted from the magnitude of the dividend (or remainder). Smaller and 
smaller increments are used until l(x 2’ ) x divisor cannot be extracted, then 

the quotient and dividend are shifted, and the system is run again to give a se¬ 
cond quotient register. The results of this subroutine are accurate to the last 
digit. 


143 THE VIBRATIONAL FREQUENCY SPECTRUM OF A COPPER CRYSTAL 

The vibrational frequency spectrum for a-iron has been calculated from the secular 
equation given below, using the atomic-force constants determined by H. Curien. This 
calculation was performed for 1784 wave vectors uniformly distributed throughout 1/48 of' 
the Brillouin zone. The WWI arithmetic program developed previously for the spectrum 
calculation of copper (face-centered cubic - F.C.C.) was appropriately modified to yield 
solutions for the body-ccntcred cubic (B.C.C.) secular determinant of d-iron. 

2 


'11 


<q> * L 


| D(q) - l(JJ | = 0 
Za( 1 - cos Uj cos u^ cos u 3 ) ♦ a 


h/3 (sin 2 u^ *■ sin^Uj) 


♦ a "( 1 - cos 2u 2 cos 2up *^"(2 - cos 2u ( cos 2u^ - cos 2UjCos 2u^) 
D^(q) - jsin u^ sin Uj( 2^3 cos u j t 4 y" cos u^ cos u^) 

Uj = 2tTaq j, u^ = ZTTaq^ . = 2iraq^ 


Curien. H. , Thesis, L'Universit^ de Paris, 1952. 
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| q I * ^ q s »qi ♦ jq 2 ♦ k qj 

1st neighbor: o * 1. 17 x 10 4 dynes/cm - 1. 19 x 10 4 dynes/cm 

Znd neighbor: 0 ' s 1. 16 x 10 4 dynes/cm s -0.25 x 10 4 dynes/cm 

3rd neighbor: £ = 0.08 x 10 4 dynes/cm Q " = y " = 0. 26 x 10 4 dynes/cm 

The composite spectra for the F.C.C. and B. C. C. cases exhibit the same general 
features. However, it is interesting to note that the two "transverse" branches for the 
B. C. C. case of a -iron show a somewhat Gaussian-like distribution, whereas for the F.C.C. 
case of copper these branches have tails only at the low-frequency end, cutting off sharply 
at the high-frequency end. One suspects that such behavior is caused primarily by the crys¬ 
tal-lattice symmetry rather than by the particular values of atomic-force constants. 

This present calculation for a body-centered cubic and the latter calculation for a 
face-centered cubic* contain no assumptions regarding the atomic forces. The force con¬ 
stants are tensor quantities and we re expe rimentally determined from thermal diffuse X-ray 
scattering experiments. All previous spectrum calculations assumed central forces and 
ron«idered only first and second nearest neighbor interaction. However, the vibrational 
spectra calculated from such a model, which is a rather special case, show no essentially 
different features from the general cases treated by this investigator. 

This work was carried out by E. H. Jacobsen of the MIT Physics Department. 

147 ENERGY BANDS IN CRYSTALS 

Dr. D. Howarth of the MIT Solid State and Molecular Theory Group has calculated 
the lower energy levels of copper for all typical symmetry points in reciprocal space; he 
has done this three times, once for each of three different one-electron potentials for cop¬ 
per. The results of these calculations are still being examined, so detailed conclusions 
cannot yet be drawn. However, the results do show that the energy value of a level depends 
more strongly than one might suppose upon the choice of the potential. It has also been 
shown that with the method used calculation of the energy at points of no symmetry in re¬ 
ciprocal space is possible although not practical; such a calculation requires over an hour 
of computer time. 
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A somewhat surprising result was the disagreement of the present method with Dr. 
Howarth's previous application of the cellular method to copper. * It was found that the rel¬ 
ative positions of the twofold and threefold degene rate levels of "d"-like symmetry were 
interchanged in the two methods, the present method agreeing roughly with those of other 
workers. Since the present method demands a slightly different form of the potential than 
docs the cellular method, one may assume that this it the origin of the discrepancy. To 
test this assumption a program was written tocalculate these energy levels by the cellular 
scheme but using the modified potential of the present method. Unfortunately, Dr. Howarth 
had to return to England before the program had been completely tested. 

149 digital methods of detecting signal from noise 

The study of digital methods of detecting sigral in noise has been completed by Dr. 

G. P. Dinneen of Lincoln Laboratory. A small amount of computer time was used at the 
beginning of this quarter to complete the study. No significant changes in the results de¬ 
scribed in Summary Report No. 38 were made. It is interesting to note that in the one case, 
that of the success-run observer, where partial theoretical results are obtainable, the cal¬ 
culations obtained by the Monte Carlo method on Whirlwind are in excellent agreement with 
these theoretical results. 

155 synoptic climatology 

During the past quarter, progress has been made, under the direction of Professor 
T. F. Malone of the MIT Meteorology Department, along six avenues of research which ^ 
follow directly or indirectly from the statement of the problems given in Summary Report 
No. 38. Each of these six points will be discussed below. 

1. M atrix Programming 

The project has advanced to a point where the solution of large-sized matrices is 
necessary. A program that will handle asymmetrical matrix, using the Crout method, has 
been written and used extensively during the past few weeks. This is a general program 
which can readily be adapted to any specific problem. The largest matrix that the routine 
will handle is determined by the drum storage available for holding the elements of the ma- 
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tr»x. I he program has already been used lor matrices up to order 3b. and it is expected 
that it will be used on matrices of theordei of 78; the largest matrix that is presently pos¬ 
sible appears to be of order 100. 

_ Pressure Prediction on New Data 

The use of the linear operators derived for pressure prediction was discussed in 
Report No. 38. Complete pressure predictions have been made for the month of January 
1951, this month was not used for the derivation of the linear operator. Twenty-four-hour 
predictions were made for each of the 91 points on a grid covering most of North America. 
Analysis of the results of these predictions has not been completed, but a preliminary ex 
amination indicates that the accuracy is rather good. This method appears to offer an ob¬ 
jective procedure for pressure prediction on a routine basis. 

3. Study on Extended R ain Forecasting 

Considerable work was done during the past quarter in an attempt to assess the in¬ 
formation available in five-day mean circulation patterns concerning the associated five- 
day mean anomalies of temperature and the five-day accumulation of rainfall. This is es¬ 
sentially an attempt to adapt the techniques embraced in this work to methods of five-day 
forecasting, as now practiced. In addition to the association of actual weather w'ith circu¬ 
lation patterns, an effort is being made to determine the predictability of the weather ele¬ 
ments. It is expected that this work will be completed during the next quarter. 

4. Precipitation 

Considerable machine time has been used in an attempt to see if it is possible to re¬ 
place a large number of circulation parameters by a relatively smaller number of param¬ 
eters which characterize the humidity distribution. The purpose of this work is to pre ¬ 
dict areal averages of 24-hour accumulation of precipitation. The work has been proceed¬ 
ing by successively more difficult stages, and a final answer on whether the moisture dis¬ 
tribution should be introduced explicitly or implicitly, in the form of circulation param¬ 
eters, has not yet been answered. 

5. Four-Dimensional Climatological Model 

This work has been proceeding along the lines indicated in Summary Report No. 38. 
Essentially it involves putting together upper-level circulation patterns with surface patterns 
and taking account of the time variation in these patterns to predict rainfall and tempera- 
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ture. A good deal of time has been spent in preparing parameter tapes for this work, but 
no actual runs have been made. 

6. Extension of Grid Size 

It has seemed worthwhile to investigate the feasibility of extending the present grid 
which is now bounded by latitudes 30* and 60* North and longitudes 65* and 125* West. The 
necessary data to do this have been obtained from the U.S. Weather Bureau in Washington, 
and the work has now advanced far enough to expect that actual machine work will be under 
way during the next quarter. 

159 WATER USE IN A HYDROELECTRIC SYSTEM 

The bulk of the work on this problem has been completed and has been written up by 
J. D. C. Little as p^rt of a doctoral thesis in the MIT Physics Department. A few more 
calculations may be made during the preparation of papers for publication. 

The problem has been: How should the storage water of a hydroelectric system be 
used when the future river flow is uncertain "* Methods currently in use for answering the 

nii^cli nn i rp hjcp/l nn injlycpc n( PPOfrVO i? C J"* ** '*’ith J H -J « snmp4 futUT? flc.V. The 

approach used here is to characterize future flow by probability distributions and determine 
the water use which minimizes the expected cost of operation. 

The Whirlwind computations consist of three parts: the calculation of the probabili¬ 
ties from the historical record of flows, the calculation of tables of best water use, and the 
calculation o f the operation of the .,-Lciii with the historical flows as input. 

During the past quarter, the optimization program has been considerably improved, 
a new scope-output routine has been added to the historical-flows program, and the major 
production runs have been made. The scope output has been very successful. Originally, 
it was planned to bring out of the machine on magnetic tape certain results of system opera¬ 
tion and then plot representative data by hand. This soon proved impracticable, both be¬ 
cause of the amount of magnetic tape used and because of the hand labor required. There¬ 
fore, a scope routine was written which displays a set of four (actually, any number) curves 
on the scope. For a given method of determining water use and for a given year in the his¬ 
torical record, curves were displayed to show the historical flow, recommended storage- 
water use, required supplemental generation, and reservoir elevation, all as functions of 
the time through the year. Although too gross to show up fine differences of operation, the 
curves give an excellent feeling for how the system is behaving macroscopically. 

Operation as determined by the expected-value method has been compared to a more 
c cnvcn*t onoperation nn •» -lmple model of a hydro system idealized from the Grand Cc lee 
plant on the Columbia River. Differences between various kinds of operation were found to 
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be small. However, the expected-value method gives a sensible stable operation with an 
average cost comewhat lower than that of conventional operation. 

166 CONSTRUCTION AND TESTING OF A DELTA-WING FLUTTER MODEL 

The problem and computational procedures of designing a delta-wing flutter model 
to simulate a given setof flexibility influence coefficient of a full-scale wing were described 
in Summary Report. 37 and 38 by M. M. Chen and S. I. Gravitx of the MIT Aeroelastic and 
Structures Research Laboratory. 

The entire programming of the iterative procedure for a three-bay structure has 
been completed and is now in working order. Four basic tapes aretnvolved In the computa¬ 
tion. The first one is a Flexo lape which is designed for the direct readin of the given in¬ 
fluence coefficients and the first estimates of the stiffnesses of each individual member. 
The second tape is in binary (556) form andcontains the geometry of the system. The third 
tape is also in binary form and corresponds to the calculation of AC's. The influence coef¬ 
ficients are computed on the basis of the first estimates of the stiffnesses, and the differ¬ 
ences are printed out. The last lape, in binary form, corresponds to the formulation of 
AX's. If the trial influence coefficients are notclose enough to the given ones, corrections 
must then be maHp tn rhe rri:*! stiffnesses. 

The program works as follows: For the first cycle all four tapes are read into the 
computer. The preliminary output from each cycle is stored on the drum right over the 
corresponding output from the previous cycle. Therefore, for each cycle of the iteration, 
only the last two computational tapes need be repeated. 

The last step of the fourth tape make, use ofCrout's matrix-inversion routine which 
solves the set of linear simultaneous algebraic equations to yield the corrections. Unfor¬ 
tunately. the principal difficulty in the iterative procedure arises in this routine; unreason¬ 
able answers were obtained because of the ill-conditioned matrix of coefficients. At pres¬ 
ent satisfactory results for a three-bay model were obtained by a trial and error scheme 
The total deviation of the flexibility influence coefficients from the given ones is on the av¬ 
erage about 7%. and the error of the influence coefficients of the diagonal element runs from 
0. 1% to 5.8%. 

During this period, a three-bay. 5/32-inch aluminum model was built and tested 
The experimental results agreed fairly well with the theoretical influence coefficients, and 
the average error is about 3%. 


167 TRANSIENT EFFECTS IN DISTILLATION 


Work has been continued by J. F. O'Donnell nf the MTT Chemical Fngineerine De¬ 
partment on two basic problems which have been described in preceding reports. The re¬ 
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suits will form a part of his Sc.D. thesis, which is bring done under the »upc rviaion of Pro¬ 
fessor E. R. Gilliland. Each of these problems involves transient effects of column hold¬ 
up in binary batch distillation. The problems are batch distillation, product takcoll at 
constant reflux ratio, and transients in continuous distillation. 

One new program has been developed and operated successfully on Whirlwind dur¬ 
ing the past quarter. It solves the equation for the limiting cas % of product takeoff in bin¬ 
ary batch distillation at constant reflux ratio, assuming the holdup in the column is neglig¬ 
ibly small. 


The basic relationship for this case is. 




( 1 ) 


where W is the fraction of the charge not yet distilled, is the distillate composition, and 
is the still-pot composition From a combination of material and energy balances and 
equilibrium relationships, can be expressed explicitly in terms of X^. However, the 

only way can be obtained from X^ is by iteration. The initial value of X^ correspond¬ 
ing Lu ^ is calculated by iicmiiuu. i'iie i callc i , as uiuscii auu X^ calculated. inr 

resulting increments in X^ are unevenly spaced. One of the problems was the question of 
how to carry out the integration. 

Initially it was hoped that by using the trapezoidal rule for the integration, solutions 
of sufficient accuracy could be obtained with reasonable machine time. Two sample cases 
were run on ihe machine, each with at least three different increment sizes. Then it was 
determined that a plot of the result, W. for any value of X w against h^. where h is the in- 

w i 

cremcnt size, was a straight line. This implied that Richardson's h extrapolation might 
be used. 

Briefly, the following procedure is used for the integration. Starting at one point, 
values of W arc calculated using two different increment sizes alternately with one equal¬ 
ing half the other. That is, two steps arc taken using the s.mallrr increment size and then one 
step with the larger. Thus, at that new point there arc available, at a single value of X^, 
values of W calculated using the two increment sizes. Then the mathematical equivalent of 
Richardson's h^ extrapolation is done, giving a corrected value of W. This is printed out. 

The results from this procedure were very satisfactory. Less machine time and 
output were used than would have been necessary if the two calculations had been done sep¬ 
arately (because it was not necessary to calculate X^ from X^ again fo» the larger step). 
Also less machine time was used than would have been necessary if each case had been done 
onlv one increment size, small enough to insure sufficient accuracy. A check was 
made each lime the extrapolation was carried out to see tha* it was not too large. 
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Using this program a series of cases has been solved, giving the product composi¬ 
tion as a function of the fraction of the charge distilled. These cases with holdup neglig¬ 
ible include all the sets of parameters which had been used previously to obtain solutions 
for the same ope ration as suming vary ing amounts of holdup present. Comparison of the two 
types of results has been valuable in gaining an understanding of the etfects of holdup. 

The major part of the work done on transients in continuous distillation has been on 
the problems of quantitative correlation of the data tor the case of a step change in feed 
composition. Two correlations were developed for the data obtained previously. Then ad¬ 
ditional solutions were made using Whirlwind to determine the effectiveness of the correla¬ 
tion* in predicting results for different parameter combinations. The work showed that the 
correlations were quite successful fo* conditions within the range previously studied. How¬ 
ever. some extreme sets of parameters gave results that were not predicted by the cor¬ 
relation satisfactorily for engineering purposes. J. 1. Iordan has submitted a Master's 
thesis to the Chemical Engineering Department on this phase of the subject. * A consider¬ 
able improvement in understanding of the situation has been achieved. 

Work is continuing on the over-all project. Though some additional data will be 
taken, the major effort will be on the utilization of available data to reach qualitative and 
quantitative conclusions. 

169 UTILIZING A GENERAL-PURPOSE DIGITAL COMPUTER IN 

SWITCHING-CIRCUIT DESIGN 

As this is the final report on this problem, a brief review of the problem and the 
results that have been obtained will be given. 

The problem undertaken by E. C. Hoy of the MIT Electrical Engineering Depart¬ 
ment involves the simplification of Boolean functions or switching functions, which are Bool¬ 
ean functions interpreted as representations of scries-parallel switching circuits. Because 
Boolean functions may be interpreted as switching-circuit representation, the problem of 
their simplification is of more than academic interest. The arbitrary criterion for sim¬ 
plicity is the fewness of independent variable appearances in a function form.* In general 
for a given function there is no unique form of any chosen degree of simplicity. Thus there 
may be several forms of the greatest degree of simplicity. It is desirable to locate at least 
one such simplest form with full cognizance of the fact that such a form has been obtained. 

A program for simplifying switching functions has been developed based on the ap- 

1. Jordan, J. I., Jr.. "Transient Conditions in Continuous Distillation Columns, " S.M. 
Thesis, Chemical Engineering Department, MIT (1954). 

1. Other criteria may be chosen, e.g., fewness of the number of variable appearances 
plus terms. 
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plication of several logical (or Boolean) identities In the course of the programming, meth¬ 
ods for adapting the digital computer to Boolean algebra were developed, including a meth¬ 
od for representing the function in storage and methods for performing the three major log¬ 
ical operations. * Input and output to the program are made in notations familiar to those 
working in the » witching-circuit field. Delayed output by means of magnetic tape is used to 
save computer time. The program requires 506 registers of fast storage in addition to the 
storage required for the statement of the particular problem. At present the program can 
effectively handle any function of nine variables. A function of any number of variables up 
to and including 15 can be handled providing its statement does not exceed the storage avail¬ 
able. By using auxiliarv storage or by sacrificing some of the input and output convenien¬ 
ces, the program can he revised to handle more functions than it does at present Provid¬ 
ing the number of variables does not exceed 15, the number of variables docs not directly 
affect the operating time of the program. Starting with standard normal forms, ^ the pro¬ 
gram required less than L i seconds to simplify 26 functions averaging over 10 terms each. 
However, although the computer -implifies functions much faster \aoout 100 times faster) 
and more accurately than people using the same simplification method, the simplification 
method that has been programmed docs not insure simplest solutions and docs not indicate 
the nature ot the solution. 

Various methods for modifying the existing simplification method were studied, and 
different approaches to the problem were sought. The most promising of the methods in¬ 
vestigated is that proposed by W. V. Quine (see reference below). His method has a proof 
indicating the nature of the solution. 

The answer to the question of whether it is feasible touse a general-purpose digital 
computer to simplify switching functions is necessarily a subjective one, depending upon 
the particular engineering problem of the designer and the equipment available to him at 
reasonable cost. 

For a more complete discussion of this work than presented here, one may refer to 
a forthcoming Electrical Engineering S.M. thesis by E. C. Hoy entitled, "Simplifying Switch¬ 
ing Functions Using a General Purpose Digital Computer." 

References: Shannon. C. E. . "A Symbolic Analysis of Relay and Switching Circuits,’* Trans¬ 
actions of the AlEE, Vol 57. pp 13-23, 1938 

Quine, W. V., "The Problem of Simplifying Truth Functions, " American Math¬ 
ematical Monthly, Vol 59, pp 521-53 1, 1952 


1. See Quarterly Progress Report No. 11 of the Project on Machine Methods of Computa¬ 
tion, p. 9. 

2. I.e. , each variable appears once in each term. 
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172 OVERLAP INTEGRALS OF MOLECULAR AND CRYSTAL PHYSICS 

The routines developed by F. J. Corbato of the MIT Physics Department and des¬ 
cribed in Summary Report No. 38 have been used lo obtain the integrals arising in the 
tight-binding calculation of graphite. These integral values arc the input data for a com¬ 
bined generation/secular-equation program which for a series of reciprocal lattice vectors, 
k, generates the dependent matrix elements of the equation, 


I H X = X S. X. E 
j ij jm ^ ij jm mm 

and solves for the eigenvalues, E mm » and eigenvectors, X^. This program was used 
for several sets of possible integral values. The results revealed that the secuiar-equation 
solutions were sensitive to inaccuracies in the S.^ matrix and that a greater number of in¬ 
tegrals than had been suspected heretofore were required for the necessary accuracy in 
the Sij matrix elements. Inasmuch as the present generation/secular-equation program is 
not easily extended to a greater number of initial integral values, the entire program is 
currently being rewritten. 

During the latter part of this quarter, time was spent consolidating and partially 
generalizing the integral-evaluation program. In addition, on the basin nf the evnerience 
gained from the present work, several improved subroutines were offered to the subroutine 
library. 


173 COURSE 6. 537 -- DIGITAL COMPUTER APPLICATION PRACTICE 

Two of the students who were enrolled in the MIT Electrical Engineering course 
6. 537, entitled "Digital Computer Application Practice," given by Professor C. W. Adams, 
have continued to work on special problems. 

The time spent on Whirlwind by H. C. Kreide has been used to develop a date- 
demonstration routine as a special problem for course 6.68, The purpose of the demon¬ 
stration routine is threefold: (1) to supply the day of the week that specific dates fall on; 

( 2 ) to supply all of the dates in a particular month upon which a given weekday falls; and 

(3) to supply the date of Easter for any desired year. 

The routine is based on the Gregorian Reform and handles all dates between 1600 
and 16,000 A. D. Calculations are based on the theory of congruences. Basis for correct 
dates and weekdays has been Cunningham's "250 Year Calendar." Basis for correct dates 
of Easter has been the 19 54 World Almanac. All tests made on the final routine have yielded 
correct results (only the years between 1753 and 2100 are given in publications). 

As it presently stands, the routine consists of four tapes, the last of which contains 
the dates about which information is requested. It will correctly answer questions about 
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dates when given in proper form. A few revisions will be made to allow the information 
given on Flexowriter tape to be more flexible and to reduce considerably the possibility of 
misinterpreting it. Also, some revisions will be made in information that is printed out 
when illegal or nonexisting dates are given the computer. Any furthei revisions Arill be to 
eliminate redundancies. Even without revisions, the program is in good operating order. 


The problem considered by S. Pet rick was the adaptation of a form of the game 
"battleship" to regular WWI code. The user of the program decides on how he wishes to 
place his five "ships" on a 16 x 16 grid. The ships have lengths from two to six squares 
long and must lie either horizontally or vertically. In addition, no two ships can be adja¬ 
cent along a line. 

1 he machine picks a square and then consults a stored table to see if a hit. has been 
made. If so, control is transferred to a polishing-off routine. The method used is the 
following; Consider square x.^. Determine the number of different ways in which a class- 
two ship can be placed on the board at this time. Repeat this process for class-three, 
-four, -five, and -six ships. Add the six resulting quotients to give an "index of shoota- 
bility" for square x^. Repeat this process for all other squares, and then pick the square 
whose index »j greatest. 

At present the routine has been found to work except for the action of the polishing- 
off routine on a class-six ship. It is hoped that this error can be found and that no other 
shortcomings ot the polishing-ofl routine will turn up. 

Only one case which worked satisfactorily has been tried, and the same distribution 
of ships was given to several human players. While Whirlwind did not take fewer guesses 
than the best human guessers (two of IS had lower scores), its score was much below the 
average human score. 


174 TIGHT-BINDING CALCULATIONS IN CRYSTALS 

The unperturbed wave function for a crystal is approximated by a linear combination 
of atomic functions. A Hamiltonian matrix is set up between these functions; the eigenvalues 
of this matrix are the energy levels of the electrons. The diagonalization of the Hamilton¬ 
ian has been done on WWI. 

In particular, this program has now been carried out for the nickel metal. This 
involved the diagonalization of 680 five by five matrices. The matrix elements were cal¬ 
culated in the machine instead of just reading in the matrix elements after computing them 
by hand. This calculation has now been completed, and the results have been processed to 
supply a density of states for nickel. This density of states represents the number of 
energy levels in the nickel metal which are in a given energy range and is important in 
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studies of magnetism and specific heat in solids. The results are now being prepared for 
publication in The Physical Review by Dr. G. Roster of the MIT Solid State and Molecular 
Theory Group. 

175 IMPURITY LEVELS IN CRYSTALS 

The calculation of impurity levels in crystals involves the solution of difference 
equations which are solved by a Green's function method. This involves finding the inverse 
of matrices (H-E) where H is the unperturbed Hamiltonian anti E is the energy of the im¬ 
purity level. The inversion is being done on Whirlwind I by Dr. G. Roster of the MIT Solid 
State and Molecular Theory Group. 

The routine lor inverting the difference of the Hamiltonian matrix and the energy E 
has been applied to chromium in order to study the energy levels of impurities in this metal. 
For each Hamiltonian matrix H, 16 different E values were subtracted from it before in¬ 
version. The total number of Hamiltonian matrices which were considered was .‘.JO. (The 
method of determining the energy levels of the impurity atom in the solid from the inverses 
of (H-E) is described by G. F. Roster and J. C. Slater in Phys. Rev. 95. 1167 1954).) 

A knowledge of the impurity levels in chromium is useful in studying thr .! nation 
of the saturation magnetic moment in the transition elements as it varies with impurity 
calculation. The work on this problem, which was to be carried out on Whirlwind, and the 
results of the computation are now being processed by the method described in the refer¬ 
ence above to give the actual levels in the solid. 
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TRANSIENT TEMPERATURE OF A BOX-TYPE BEAM 


The time histories of temperature and stress at several points along the periphery 
of a thin-walled rectangular cross section have been calculated on WWI by L. Schmit of the 
MIT Aero-Elastic and Structures Research Laboratory. By using an extruded rectanoular 
cross section (of 63X-T5 aluminum alloy), the joint contact admittance problem has been 
avoided. The heat input is constant with time but varies chordwise. All spanwisc effects 
are neglected. The temperature gradient across the thin wall thickness is assumed to be 
negligible. The specific heat, thermal conductivity, and weight density are assumed not 
to vary with temperature. However, the modulus of elasticity and the thermal coefficient 
of expansion are assumed to vary with temperature. The temperature distribution is 
carried out by an uncoupled finite-difference procedure. The physical grid employed con¬ 
sists of 54 1/8-inch square elements. A time increment of 0.02 second is selected (M e 
q —-- = 6.25). Upon completion of every 50 time cycles the stress distribution is deter¬ 
mined on a basis of the current temperature distribution. The stress calculations are 
based on an extension of the elementary-beam theory to include thermal effects. Both the 


28 


C RELEASE. CASE 06-1104. 


MATHEMATICS. CODING. AND APPLICATIONS 


temperature and time histories obtained for a heating duration of 25 seconds are to be 
compared with experimental results. 

180 CROSSCORRELATION OF BLAST-FURNACE INPUT-OUTPUT DATA 

A combined program which will compute and plot for any given functions a corre¬ 
lation function and its Fourier transform is now being optimized and tested. The program 
is designed to minimize the tape-prepa ration time, operating complexity, and total com¬ 
puter time involved in correlating a pair of functions and computing the transform of the 
cor relation function. In addition, the program eliminates all intermediate hand computation 
and tape preparation, since a single continuous run completes all computation. 

A further aid in operating the program is the use of only the oscilloscope camera 
as the output device; the problem of efficiently scheduling magnetic-tape utilization in run¬ 
ning a series of problems is greatly simplified. 

The program will handle lunctions up to 250 values in length and uses the Fourier 
transform program developed under Problem 171 by D. T. Ross. This program has the 
feature of automatically reducing spurious Gibbs ripples in the spectra. Both cosine and 
sine transforms are computed. However, when a function is autocorrelated, only the co¬ 
sine transform is required, and the program may be stopped after this result is obtained, 
thus eliminating a lengthy and unnecessary computation. 

It is hoped that the availability of such a program in the Laboratory files will ma¬ 
terially reduce the programming effort involved in many types of problems, since the 
processing which this program automatically carries out has very wide applications in 
scientific and engineering computation. This work was initiated by R. G. Mills for an S. B. 
thesis in the MIT Electrical Engineering Department. 

183 BLAST RESPONSE OF AIRCRAFT 

During this past quarter, three new phases of the blast loading problem have been 
programmed, and test runs have been made on WWI. 

The problem, being programmed by Y. Shulmanof the MIT Aero-Elastic andStruc- 
turcs Research Laboratory, is to determine the dynamic response of an idealized airplane 
upon encountering a symmetric blast gust. The airplane is approximated by one fuselage 
mass and twe wing masses connected by two uniform massless beams. The method of 
solution is, in general, similar to that of the cantilever-beam problem reported earlier. 

Unlike the cantilever-beam case studied earlier, the structural characteristics of 
the uniform beam, which represents the wing structure, are represented by a moment- 
vs. -curvature curve having a discontinuity in bending moment at the instant that *.vir.g 
buckles. Sub* , **q , i*»r'* buckling, the ;tf istoiice oi the beam decreases as its root curvature 
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increases. In view of this consideration, the programming of the problem provides an 
additional feature: namely, the maximum response of the structure is defined either as 
the actual peak curvature of the curvature at which the resistance of the beam is zero. 

The blast force acting on the wing is formulated according to unsteady aerodynamic 
theor r .ind lumped at the mass stations. The exponential form of Wagner's function is 
used, and hence the subroutine for e* is used in the program. 

In addition to the blast force and inertial force, the airplane is also subjected be¬ 
cause of motion to damping forces which were neglected in the cantilever-beam case. For 
one phase of the present problem, the damping force is formulated according to unsteady 
aerodynamic theory and reduced to a recurrent form. For a second phase, the damping 
force is assumed quasi-steady, and the numerical-integration procedure is greatly simpli¬ 
fied. 

Tapes of the above two phases of the problem have been tested and checked with 
desk solutions. Results of the runs also bhow that the assumption of quasi-steady damping, 
though satisfactory in the linear range of the structural characteristics, yields large error 
in the post-buckling range. 

Furthermore, the test runs of the tape also indicate that the machine time needed 
for the automatic reduction of the time increment is longer than desirable. Consequently, 
a convergence test has been made for a set of parameters of practical interest without 
automatic reduction. It was found that for a nondimensional time increment of 0. 1 the 
error in amplitude is within 2%, and the error for the time when the peak occurs is neg¬ 
ligible. 

A third phase of the problem is to determine the peak blast-gust velocity which will 
cause a given wing structure to fail at various altitudes and airplane velocities. For a set 
of parametric values, the problem thus involves a cut-and-try procedure so that the peak 
response is equal to the given failure curvature. A tape has been prepared with the pro- 
visionof an automatic increase in the value of the peak blast gust. The tape is, at present, 
being tested on WWI. 

The production runs scheduled for the near future consist of 50 runs for the first 
tape and 40 runs for the third tape. The machine time needed for the first tape is esti¬ 
mated at about 3 minutes per run. No estimate can be made at this time for the third 
tape. 

184 SCATTERING OF ELECTRONS FROM HYDROGEN 

The calculation of the two-dimen'*onal integrals described in Summary Report No. 
38 has been largely completed. Each integral corresponds to a probability of a virtual 
transition to an excited state of the hydrogen atom. These expressions were evaluated for 
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three values of the energy of the incident electron and five values of the scattering angle. 
It was possible to perform the integrations analytically when the scattering argle was zero 
degree. This exact calculation was a rheck on the numerical procedure, and the two re¬ 
sults checked to the desired accuracy (three significant figures). 

in the case of the virtual transitions to the continuum states, it had been planned to 
reduce what we re originally three-dimensional integral s to a two-dimensional form and then 
perform the numerical quadrature, as in the case of the transitions to the discrete states. 
However, it finally had to be conceded that the approximations that led to this reduction 
were unsatisfactory, and a program was constructed to perform the original three-dimen¬ 
sional integrals numerically. These terms are also being checked against the results of 
simpler calculations for zero scattering angle. 

This work is being carried out by M. C. Newstein of the MIT Physics Department 
and will be included in his Ph. D. thesis. 

185 A SCALE OF TURbULENCE 

The correlation of eddy velocities or fluctuations between two points within a turb¬ 
ulent flow will generally vary inversely with the magnitude of the distance between the 
points. This offers a method by which eddy size may be defined, the length of the eddy be¬ 
ing that separation of the points at which the correlation becomes insignificant. 

Atmospheric turbulence covers a wide range of motions with associated length scales 
varying from large disturbances containing most of the energy of the motion to those smal¬ 
lest eddies which dissipate energy as heat by action against viscosity. J. W. Tukey states 
that the minimum resolution of fluctuations is fixed by the sampling interval on a stationary 
stochastic process. An upper limit is set by the length of the sample. 

Tukey's analytic filter was applied to seven sets of wind-velocity data. These data 
were simultaneous observations of wind speed at a height of about 2.5 meters made by four 
standard Hastings hot thermocouple anemometers with a time constant of 0.25 second. The 
anemometers were aligned either parallel or normal to the mean wind flow with a separa¬ 
tion of 1, 2, and 4 meters between individual units on five of the observations and 3, 6, and 
12 meters on the remainder. The sampling periods were from 392 to 510 seconds long. 
The sampling interval was 1 second. 

The analytic filter operates on the sequence 

*1‘ *2' " ' x i - I' x i + 1. X N - r *N 

by forming the following differences: 
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.4. Fir^t-ordcr differences of the j-th probe. 


i* 1. 2, . . 


b. Second-order differences for the j-th probe, 

d - i (x t ♦ x. ,) - (x.. , ♦ x. ) 

j j ' 4i 4i-1 ' 4i-2 4i- 3' 

c. Third-order differences for the j-th probe. 


i = 1, 2, . 


d j Ul = |< x 8l + •" + x 8i-i> - (x 8,-4 

chore represents the integral part of N/ 2^ 


*8i-7> }• iBl ’ 2 .[ N 8 ] 


According to Tukey, the sequences of differences as formed above provide scries 
in which maximum spectral power is concentrated near periods of 2^ seconds per cycle, 
where p is the order of the difference. The filter is only a rough one since it lacks sharp 
cutoff characteristics. In this problem, only the first three differences are considered. 

Other investigators have shown that the sums of squares of the differences defined 
above are proportional to the variance and energy at their associated frequencies. 

Seven parameter tapes were prepared for presentation to Whirlwind using the ofc 
served wind-speed reading. A program was then prepared in which the differences pre¬ 
viously defined were formed intermediately without printout. For each difference, for each 
parameter tape. Whirlwind was called upon to provide the following numbers as its contri¬ 
bution to the solution of the problem 

[n/2 P ] [n, 2 P ] [n, 2 P ] [n/2 P ] 

2 V I Vik- 2 W 2 -Ho 2 ** 2 - 3. « 


where d^ is the i-th difference in the sequence of differences for the j-th probe. Thus 630 
numbers were printed by Whirlwind. 

With the results from Whirlwind, certain statistical tests were made on the data in 
order to ascertain the degree of validity of the assumption that the d.^ are normally dis¬ 
tributed with mean of zero and with homogeneous variances for synoptic observations. 
This was accomplished by considering each difference for each of the seven sampling 
periods separately as a subject for the following battery of tests. First, the homogeneity 
of variances was tested by using the test described by Hocl for this purpose. Then an 
Analysis of Variance was performed to test the hypothesis that the means from each of the 
four probes were essentially the same, and the Student's "t" Test was applied to show that 
this mean was not significantly different from zero. Finally, a Chi-square Goodness of Fit 
Test was applied to the d.^ to determine if they could considered a**, having a -.ormal 
distribution with . w «*,».! variance equal to the arithmetic mean of the variances 
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of the four individual probes. 

For each of the differences for each sampling period, the correlation coefficients 
for the six possible combinations of the probes were computed. These coefficients were 
plotted against distance, and the indicated length at which the coefficient becomes sensibly 
zero was noted. Thus, an eddy with a predetermined frequency will have a length, or 
scale of turbulence, associated with it. 

In addition to the usual cor relation coefficient, the d^ were first squared (by Whirl¬ 
wind) and then a figure analogous to the correlation coefficient computed and plotted against 
probe separation. Also, the quantities I i d^ - d.^ I were plotted against probe separa¬ 
tion. The last two plotted figures will be compared to the plotted correlation coefficients 
to determine the degree of agreement between them. 

Production runs have been satisfactorily completed on all seven sets of data. The 
results have been described in an S.M. thesis written jointly by J. Howcroft and J Smith 
for the MIT Meteorology Department. 


186 TRACKING RESPONSE CHARACTERISTICS OF THE HUMAN OPERATOR 

A series of power-density and cross-power-density spectra has been computed for 
J. I. Elkind of the Lincoln Laboratory at MIT by WW1. These spectra were computed by 
taking the Fourier transform of the auto- and crosscorrelation functions of stimulus and 
response. The spectra that have been computed were obtained from an experiment designed 
to measure (1) the variability in the response characteristics of particular subjects, and 
(2) the variability in response characteristics among several subjects 

Only a few of the spectra have been used to obtain human operator transfe r functions 
and noise power-density spectra. Most of the results already computed by WWI have not 
yet been reduced to this form. The transfer functions that have been obtained show that 
the human operator acts like a low-pass filter of band width about as great as the band width 
of the stimulus he is tracking. 


190 ZEEMAN AND STARK EFFECT IN POSITRON1UM 

Positronium, the transient electron-positron atom, annihilates from S-states, with 
the transfer of the rest mass energy to the free radiation field. From the *S. states the 

l * 

annihilation process yields 3 quanta, from Sq states two. The gross energy structure of 
the atom is the same as that of hydrogen save for a constant factor of 1/2 arising from the 
reduced mass of the "nucleus. ' The fine and hyperfine structure is entirely different. 

In then = 2 group of states the fine structure pattern perturbed by magnetic and elec - 
trie fields has been computed by diagonalizing the Hamiltonian for the system. The elec¬ 
tric fields arise experimentally from the motion of the atom perpendicular to the rnspnetic 
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field in which it is formed. 1 he Hamiltonian matrix is 14 x 14 and has been diagonalized 
for values of the magnetic field of 500, 100C. 2000, 3000, 4000, and 5000 gausses and for 
the values of the motional energy of 0, 0.025, 0. z5, and 0.5 electron-volt. 

The states of n = 2 all have finite lifetimes (2- or 3-quantum annihilation or optical 
emission), and an approximation to the reciprocal lifetimes of the perturbed states F y° ta * 

(i 3 1.14) may be found from the unperturbed reciprocal lifetimes 3 1.14) 

by computing 


F Total 

* s 


£ (a ik>v 


where a^ is the perturbed eigenvector. 

The branching ratio for 2- to 3-quantum annihilation of the atom may be found by 
forming the partial sums 


singlet only 


r r- s <*u>Vj 

triplet only 


For each energy eigenvalue computed, these partial reciprocal lifetimes and the 
total reciprocal lifetime were also found. 

The present program of computations is now completed. The study was carried out 
by H. W. Kendall of the MIT Radioactivity Group. The routines for diagonalizing the Ham¬ 
iltonian and for computing the partial and total reciprocal lifetimes were developed by Dr. 
A. Meckler of the MIT Solid State and Molecular Theory Group. 


192 FREQUENCY AND PHASE SPECTRUM ANALYSIS OF SEISMOGRAMS 

The following constitutes a brief resume of the work accomplished under Problem 
192 by William P. Walsh. A more detailed account of the work may be found in the author's 
Ph. D. thesis submitted to the Department of Geology and Geophysics, 16 August 1954. 

Throughout the existence of seismology surprisingly little has been done to ascertain 
by means of exact analysis the frequency content of recorded seismic motion, in compari¬ 
son to the theoretical work which has taken place in this regard. The U.S. Coast It Geo¬ 
detic Survey has, we understand, revently carried out a good deal of work in this direction. 
Their attention was confined to strong-motion earthquakes, however. Many investigators 
in the petroleum industry have recently become interested in the frequency analysis of re¬ 
producible prospecting records. As far as we know, the aforementioned is all the experi¬ 
mental research that has been projected in this direction in the fields of exploration and 
earthquake seismology. 
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In our work we have endeavored to apply the recent findings of Tukcy (outlined in 
our thesis and initial report) to the spectrum analysis of seismograms and lo exploit as 
much as possible MIT's WWl computer to that end. 

Our purpose was not only to devise a means of seismogram spectrum analysis but 
also to apply such analysis in several cases in order to obtain useful information concern¬ 
ing surface-wave dispersion, microseisms, and direct, refracted, and reflected energy of 
earthquakes and explosion .. 

The procedure we devised consisted of computing spectra (k la Tukey) from hun¬ 
dreds of successive overlapping intervals and displaying them in the fashion of a contour 
map - - called here traveling spectra. 

Application wi the aforementioned process to several earthquake records obtained 
from the observatory in Weston, Massachusetts, revealed some interesting dispersion 
curves for both Rayleigh and Love waves over Atlantic and continental paths. Those curves 
for surface waves which traveled the Atlantic path were strikingly similar, whereas those 
curves for surface waves traveling over the continent were most dissimilar. This probably 
indicated a fair degree of lateral inhomogeneity in the continental portion of the earth’s 
crust by way of comparison to the oceanic portion. It appears that exact comparison of our 
observed curves with models which more clearly delineate the complexities of the earth's 
crust must await the calculation of theoretical curves obtained from at least a five-layered 
model (because of difficulties in these calculations only curves of two-layered models have 
thus far been determined). 

Similar techniques were also applied to some microseisms (seismograms recorded 
at Weston) due to a storm which existed over the Atlantic near the coast of New England. 
Our analysis revealed the existence of a resonance phenomenon (i.e., the depth of water un¬ 
der the meteorological "low” or cold front dictates the band width and prominent frequen¬ 
cies of the microseism spectrum) and further corroborated the theory of the origin of mi- 
croscisms proposed by Longvct-Higgins. 

Since the analysis of earthquake seismograms suggested the use of traveling spec¬ 
tra as a means of determining the travel time of reflections, several trials were made on 
three prospecting records where reflections were both visible and invisible. The results 
were rather encouraging, for all types of reflections were put into evidence quite well. 

193 EIGENVALUE PROBLEM FOR PROPAGATION OF ELECTROMAGNETIC WAVES 

The field radiated from a Hertzian dipole to points well beyond the horizon over a 
perfectly reflecting earth through an atmosphere whose index of refraction decreases lin¬ 
early with height (constant gradients) can be conveniently calculated as an eigenvalue prob¬ 
lem. The eigenvalues of the allowed modes are determined by the boundary conditions at 
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the earth's surface, at the air-vacuum interface at height g. and at heights above the air. 
The real parts of the eigenvalues determine the phase velocities of the allowed modes, ard 
the imaginary parts give the attenuation rate with horizontal distance around the earth. The 
earth radiuB and wavelengths considered, of the order of a few meters or centimeters, en¬ 
ter also into the numerical values of s and g, the index of refraction gradient and air-layer 
thickness, respectively. 1 he first step, finding the eigenvalues w of the propagation modes 
consistent with the boundary conditions, leads to the problem of finding the complrx roots 
of the following equation: * 


h|(w - . K ) 8s l iv/s 3/i 


4 3/2 

— w 
3 


where w is a complex quantity, and s and g are numerical, real constants. The functions 
hj and h^ arc described and tabulated in the book, "Tables of the Modified Hankel Functions 
of Order One-third and of their Derivatives," published by the Harvard University Pres*. 
Cambridge, Massachusetts, in 1945. 

The values of w are uniformly separated, and approximate values are knewn or can 
be obtained by extrapolation from those already computed. Starting with an approximate 
value of w, an accurate value in the range so far investigated can be obtained by two ap¬ 
plications of a computing procedure. 

For each accurate value of w there is to be found the value of the following product 
of factors: 

t 3/2 x ./2 u -»/3 ( 1 - . J 

16s*w 2 


hj(w - sg + sZj) hj(w - sg + sZ 2 ) 
h 2 (w - sg + sZ j) h j(w - sg) 


h (w - sg + sZ. 


h 2 (w - sg + sZ 2 ) hj(w - sg) 


expon iXsw-iZw' + i . 


1. For the derivation of (1) and (2) see Lincoln Laboratory Technical Report No. 38 (T. J. 
Carroll, R. M. Ring. February 12, 1954). 
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where X is the distance of transmission multiplied by a known constant factor, and Zj and 
Z^ are the heights of the transmitting and receiving antennas, respectively, multiplied by 
a known constant factor. 

It is possible toiind all the values ol w for which Expression (2) has an appreciable 
size. The corresponding values of (2) are added together as complex quantities, the sum 
being the signal strength. A curve of signal strength as a function of X. or distance, may 
be plotted, and, if desired, a curve of signal strength as a function of Z^, or height of re¬ 
ceiving antenna, may be plotted. 

Numerical computation cove ring about nine-tenths of the work of calculating Expres¬ 
sion (2) has been performed on the Whirlwind digital compute r and checked by a calculation 
using a desk machine. A program tape for calculating Expression (2) has been prepared 
a^d is being tested. 

Measurements of signal strength from the literature will be compared with the com¬ 
puted values. 

This work is being carried out by Professor H. B. Dwight of the Lincoln Labora¬ 
tory at MIT. 

194 AN AUGMENTED PLANE WAVE METHOD AS APPLIED TO SODIUM 

In selecting a set of approximate wave functions of a periodic crystal potential to 
form the basis functions of the space in which the energy is to be diagonalized, one must 
first decide from which bands these functions are to be chosen and secondly from which set 
of neighbors surrounding the central cell in reciprocal space they are to be chosen. A cus¬ 
tomary and physically motivated method of choosing the latter is this: 

The cells of reciprocal space are divided into rings of increasingly distant neighbors 
of the central cell, all cells of a ring being equidistant from the origin. Each enlargement 
of the space in which the energy is to be diagonalized which does not increase the number of 
bands represented requires adding an entire ring of neighbors. 

This space can first of all be divided into invariant subspaces, each subspace being 
characterized by a particular value of wave vector. Each of these spaces cuts across rings 
and is characterized by the wave vector of the functions which form the intersection of the 
subspace with the central cell. This wave vector is referred toas the reduced wave vector 
of the set of functions which compose the subspace. For certain wave vectors the subspacc 
can in turn be reduced into invariant subspaces each of which is characterized by a ring 
and by the manner in which it transforms under the operations which leave the subspace in¬ 
variant. This set of operations commutes with the energy, of course, and is called the 
group of the wave vector of the space. It can now be arranged that in each subspacc of the 
same symmetry, but of different rings, the basis functions be chosen so that any operation 
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the earth's surface, at the air-vacuum interface at height g. and at heights above the air. 
The real parts of the eigenvalues determine the phase velocities of the allowed modes, and 
the imaginary parts give the attenuation rate with horizontal distance around the earth. The 
earth radius and wavelengths considered, of the order of a few meters or centimeters, en¬ 
ter also into the numerical values of s and g. the index of refraction gradient and air-layer 
thickness, respectively. The first step, finding the eigenvalues w of the propagation modes 
consistent with the boundary conditions, leads to the problem of finding the complex roots 
of the following equation 


h,(w - »g) 
h g(w - *g) 


4 3/2 

_ w 

3 


where w is a complex quantity, and s and g are numerical, real constants. The functions 
hj and h^ are described and tabulated in the book, "Tables of the Modified Hankel Functions 
of Order One-third and of their Derivatives," published by the Harvard University Press, 
Cambridge, Massachusetts, in 1945. 

The values of w are uniformly separated, and approximate values are known or can 
be obtained by extrapolation from those already computed. Starting with an approximate 
value of w, an accurate value in the range so far investigated can be obtained by two ap¬ 
plications of a computing procedure. 

For each accurate value of w there is to be found the value of the following product 


of factors: 


t 3/Z x 1/2 ir l/3 ( l__^ 
16s w 


•>,(» - hj(w - »g + sZp 

h 2 (w - sg + sZ ( ) Uj(w - sg) 
hj(w - sg + sZj) h“(w - sg) 

( h 2 ( w - sg + »Z 2 ) hj(w - sg) 
hj(w - sg + sZ 2 ) h 2 (w - sg) 

export f i X s^w - i ^ w 3 ' 2 + i iZL j 


For the derivation of (1) and (2) see Lincoln Laboratory Technical Report No. 38 (T. J. 
Carroll, R. M. Ring, February 12, 1954), 
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where X is the distance of transmission multiplied by a known constant factor, and Zj and 
Z^ are the heights of the transmitting and receiving antennas, respectively, multiplied by 
a known constant factor. 

It is possible tofind all the values of w for which Expression (2) has an appreciable 
size. The corresponding values of (2) are added together as complex quantities, the sum 
being the signal **»rength. A curve of signal strength as a function of X, or distance, may 
be p.otted, and, if desired, a curve of signal strength as a function of Z^, or height of re¬ 
ceiving antenna, may be plotted. 

Numerical computation covering about nine-tenths of the work of calculating Expres¬ 
sion (2) has been performed on the Whirlwind digital computer and checked by a calculation 
using a desk machine. A program tape for calculating Expression (2) has been prepared 
and is being tested. 

Measurements of signal strength from the literature will be compared with the com¬ 
puted values. 

This work is being carried out by Professor H. B. Dwight of the Lincoln Labora¬ 
tory at MIT. 

194 AN AUGMENTED PLANE WAVE METHOD AS APPLIED TO SODIUM 

In selecting a set of approximate wave functions of a periodic crystal potential to 
form the basis functions of the space in which the energy is to be diagonalized, one must 
first decide from which bands these functions are to be chosen and secondly from which set 
of neighbors surrounding the central cell in reciprocal space they are to be chosen. A cus¬ 
tomary and physically motivated method of choosing the latter is lliis: 

The cells of reciprocal space are divided into rings of increasingly distant neighbors 
of the central cell, all cells of a ring being equidistant from the origin. Each enlargement 
of the space in which the energy is to be diagonalized which does not increase the number of 
bands represented requires adding an entire ring of neighbors. 

This space can first of all be divided into invariant subspaces, each tubspau. 'cing 
characterized by a particular value of wave vector. Each of these spaces cuts across rings 
and is characterized by the wave vector of the functions which form the intersection of the 
subspace with the central cell. This wave vector is referred toas the reduced wave vector 
of the set of functions which compose the subspace. For certain wave vectors the subspace 
can in turn be reduced into invariant subspaces each of which is characterized by a ring 
and by the manner in which it transforms under the operations which leave the subspace in¬ 
variant. This set of operations commutes with the energy, of course, and is called the 
group of the wave vector of the space. It can now be arranged that in each subspace of the 
same symmetry, but of different rings, the basis functions be chosen so that any operation 
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the earth's surface, at the air-vacuum interface at height g, and at heights above the air. 
The real parts of the eigenvalues determine the phase velocities of the allowed modes, and 
the lm.igina ry parts give the attenuation rate with ho rizontal dislanc e around the earth The 
earth radius and wavelengths considered, of the order of a few meters or centimeters, en¬ 
ter also into the nunienc.il values of s and g. the index of refraction gradient and air-layer 
thickness, respectively. The first step, finding the eigenvalues w of the propagation modes 
consistent with the boundary conditions, leads to the problem of finding the complex roots 
of the following equation: 


hj(w - sg) 
'" 2 (w - sg) 



e 


-itr/3 


3/2 


expon 



3/2 


(1) 


where w is a complex quantity, and s and g are numerical, real constants. The functions 
hj and h^ are described and tabulated in the book, "Tables of the Modified Hankel Functions 
of Order One-third and of their Derivatives," published by the Harvard University Press, 
Cambridge, Massachusetts, in 1945. 

The values of w are uniformly separated, and approximate values are known or can 
be obtained by extrapolation from those already computed. Starting with an approximate 
value of w, an accurate value in the range so far investigated can be obtained by two ap¬ 
plications of a computing procedure. 

For each accurate value of w there is to be found the value of the following product 
of factors: 


^3/2 x l/a jg-1/3 ( 1 - S 


16s 


3 

rx* 


hjlw - sg 4 sZj) h ( (w - sg 4 sZ 2 ) 


h 2 (w - sg 4 sZj) 
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1. For the derivation of (1) and (Z) see Lincoln Laboratory Technical Report No. 38 (T. J. 
Carroll, R. M. Ring, February 12, 1954). 


where X is the distance of transmission multiplied by a known constant factor, and /. j and 
are the heights of the transmitting and receiving antennas, respectively, multiplied by 
a known constant factor. 

It is possible to find all the values of w for which Expression (2) has an appreciable 
j;ae. The corresponding values of (2) are added together as complex quantities, the sum 
being the signal strength. A curve of signal strength as a function of X, or distance, may 
be plotted, and. if desired, a curve of signal strength as a function ol Z^, or height of re¬ 
ceiving antenna, may be plotted. 

Numerical computation cove ring about nine - tenths of the work of calculating Expres¬ 
sion (2) has been performed on the Whirlwind digital computer and checked by a c alculation 
using a desk machine. A program tape for calculating Expression (2) has been prepared 
and is being tested. 

Measurements of signal strength from the 1 iteraturc will be compared with the com¬ 
puted values. 

Ihis work is being carried out by Professor H. B. Dwight of the Lincoln Labora¬ 
tory at MIT. 

194 AN AUGMENTED PLANE WAVE METHOD AG APPLIED TO SODIUM 

In selecting a set of approximate wave functions of a periodic crystal potential to 
form the basis functions of the space in which the energy is to be diagonalized, one must 
first decide from which bands these functions a»-e lobe chosen and secondly from which set 
of neighbors surrounding the central cell in reciprocal space they are to be chosen. A cus¬ 
tomary and physically motivated method of choosing the latter is this: 

The cells of reciprocal space are divided into rings of increasingly distant neighbors 
of the central cell, all cells of a ring being equidistant from the origin. Each enlargement 
of the space in which the energy is lobe diagonalized which docs not increase the number of 
bands represented requires adding an entire ring of neighbors. 

This space can first of all be divided into invariant subspaces, each subspace being 
characterized by a particular value of wave vector. Each of these spaces cuts across rings 
and is characterized by the wave vector of the functions which form the intersection of the 
subspace with the central cell. This wave vector is referred toas the reduced wave vector 
of the set of functions which compose the subspace. For certain wave vectors the subspace 
can in turn be reduced into invariant subspaccs each of which is characterized by a ring 
and by the manner in which it transforms under the ope ration s which leave the subspace in¬ 
variant. This set of operations commutes with the energy, of course, and is called the 
group of the wave vector of the space. It can now be arranged that in each subspace of the 
same symmetry, but of different rings, the basis functions be chosen so that ».«y np ... '.on 
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oi the group of the wave vector has the same representation in each space. When the en¬ 
ergy is diagonalized only wave functions with the same reduced wave vector and which 
transform equivalently will combine to form the new basis functions of the space. Thus if 
the »e reductions can be effected, a much smaller secular problem need be solved than if 
one disregarded the invariance of the subspaces and were content to let the secular prob¬ 
lem do the job. To effect this reduction the proper basis functions must be set up in each 
successive ring which is added to the problem. 

Mr. M. M. Saffren of the MIT Solid State and Molecular Theory Group is testing a 
program which will first determine the intersection of the space characterized by a wave 
vector with the various rings ar.d then will encode the proper coefficients with which to 
combine functions in the inter section toobtain a subspace of the proper symmetry. The pro¬ 
gram is restricted to handle crystal* only of the cubic group (this includes face- and body- 
centered crystals). The program is to be first added to Dr. Howarth's program for cubic 
crystals and will completely mechanize that program, the input becoming extremely simple. 
The routines of the program being tested will also be used in a program which is to calcu¬ 
late the energy bands in sodium according to the method of symmetrized augmented plane 

1 

waves. 

195 INTESTINAL MOTILITY 

The analysis of intestinal-motility records being carried out for Dr. J. T. Farrar 
of the Gastroenterological Section of the Evans Memorial Hospital involves the autocorre¬ 
lation and transformation of 80 sets of data. This work has been described under Problem 
107 in Summary Reports No. 37 and 38. 

Twenty-four sets of data have been autocorrelated using a total time of 46 minutes 
with 15,900 points input, 1414 points output (101 points output per run). The average run 
has 700 points of data, giving an average of approximately i minutes per run. 

These statistics indicate that, because of the high speed and efficiency of WWI, auto¬ 
cor relations obtained by these methods are not only more accurate but also compare favor¬ 
ably in cost with autocorrelations obtained from special-purpose analog computers. Addi¬ 
tional economies in obtaining power spectra will be real! Led, since the punched-tape output 
of the autocorrelation program serves as input to the Fourier transform program with¬ 
out further processing. 


196 SINGLE-ADDRESS COMPUTER 

The single-address computer (SAC) was developed for use in the MIT 1954 summer 
bt»iiontour»r 6.531, OAL. is a modification oi the SS computer, which was completed ear¬ 
lier under Problem 140. A detailed description of the SS computer was presented in Sum¬ 
mary Report No. 35, Section 5. 4. The features which distinguish SAC from the SS computer 
are the following 

Read-ln and Conversion 

A SAC program tape, punched by Mcxowriter equipment, is read in and converted 
to WWI binary information by the SAC conversion routine. A new feature which has been 
added is the option of erasing SAC storage when reading in. Erasure is effected by pres¬ 
sing the lower activate (erase) button before pressing the read-in button. This causes SAC 
quick-access storage and the simulated magnetic tapes (se: below) to be cleared before the 
program tape is converted. * Failure to press the erase button before reading in causes the 
program information to be superimposed on the previous contents of SAC storage without 
disturbing the information on the magnetic tape. 

The SAC Magnetic Tapes 

An auxiliary-storage medium, consisting of four magnetic-tape units, has been in¬ 
corporated in SAC. Each of the units has a capacity of 990 SAC words. Words are always 
transferred between SAC storage and magnetic tape in blocks of ten. Thus, each tape is 
divided into 99 blocks, numbered from 1 to 99. 

Three instructions are used in coding for the magnetic-tape units. These provide 
tape-search, tape-read, and tape-write facilities. The tape-search instruction (mts xyz) 
is used to select one of the four tape units and to position it before the desired block; the 
tape-read (mtr alb) and tape-write (mtw alb) instructions operate on the next block in se¬ 
quence on the most recently selected tape unit. 

The magnetic-tape information is stored in WWI not on the magnetic-tape units but 
on the buffer drum. This leads to faster SAC operation and facilitates the WWI coding; the 
external behavior of SAC simulates in all respects the use of true magnetic-tape equipment. 


-- 1. Provision has been made for automatically resetting the SAC tape units to contain pre- 

1. Saffren, M. M. ; Quarterly Progress Report, No. 13; July IS, losa Solid State a;,J Mo- determined initial data whenev-r this i. desired. Instead of being erased, 

lecula: Theory Group, MIT; no 44 ao. 


38 


39 







APPROVED FOR PUBLIC RELEASE. CASE 06-1104. 


MATHEMATICS. CODING, AND APPLICATIONS 


MATHEMATICS. CODING. AND APPLICATIONS 


Modification of Instructions 

The scope of the arithmetic operations in SAC has been broadened to facilitate the 
modification of instructions during the course of a program. The cycle-count lines greatly 
reduce lhe need lor changing instructions in storage, but in certain cases the ability to make 
such changes is desirable. In the SS computer, operations combining instructions and inte¬ 
gers are not permitted. In SAC, such operations may be executed, provided that the result 
is a proper SAC word. In particular, an integer may be added toor subtracted from an in¬ 
struction to produce an instruction with a modified address, and iwo instructions with the 
same operation section may be subtracted to produce an integer equal to the difference of 
the address sections. * These procedures are similar to operations normally executed by 
many digital computera (including WWI) and add considerable flexibility to the treatment of 
instructions in SAC. 

The Instruction ripO 

The SAC instruction code contains one additional instruction, ripO, not found in the 
SScodc. This instruction permits the programmer tocall for the read-in of a new program 
tape. When ripO is executed, control is transferred to the SAC conversion routine. The 
instruction is equivalent, in effect, to pressing the read-in button without pressing the er¬ 
ase butten. The tape read in using the ripO instruction must therefore be complete within 
itself; no reference may be made to symbolic addresses which were assigned in preceding 
program tapes unless these addresses are reassigned in the new tape. 

The Permanent Table of Constants 

There is available in SAC a setof fixed quick-access registers whose contents may 
be read but not changed by the routine. Register 0 contains the integer 0. and registers 
290 to 298 contain in sequence the integral powers to ten from 10® = 1 to 10** = 100,000,000. 
These constants may be utilized by the programmer i-. exactly the same manner as integers 
stored within his routine, except that any attempt to modify them will result in a post mortem. 

The Treatment of "Excess" 

Arithmetic operations on integers in SAC may yield results which are too large to 
be stored in a SAC memory location. To allow for the generation of such integers, the SAC 


1. In all cases, the counter letter, if any, is ignored. 
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accumulator is made double the length of a memory location. The part of the accumulator 
which contains any overflow that may be generated is called the "excess." Any attempt to 
transfer the contents of the accumulator to storage when the contents of excess are not zero 
will result in a post mortem. 

In the SS computer, the instruction txi air causes the contents of excess to be trans¬ 
ferred to register ale and excess to be cleared. The contents of excess, however, are not 
simply related in the decimal system to the integer which originally was generated. * Con¬ 
sequently, the use of the txi instruction is in most cases quite awkward, and it is frequent¬ 
ly more convenient to clear excess using an arithmetic division by a suitable power of ten. 
This is analogous to the decimal-shift-right found in many decimal digital computers. The 
txi instruction has been abandoned in SAC in favor of the more easily understood and more 
logical arithmetic division. The jix ale instruction, used to determine whether excess is 
occupied, has been retained. 

Timing of SAC Programs 

At the termination of a SAC program, either by the instruction StpO or by a compu¬ 
tation post mortem, a fictitious time is automatically printed by the computer, indicating 
how long a physically existent SAC would have required to execute the program performed. 
This time is computed during the course of the program, giving a preassigned weight to 
each instruction as it is carried out. The weights assigned to the instructions are based on 
1 millisecond for most operations, with proportionally greater times required for multipli¬ 
cation, division, and the operations using input and output equipment. 

In calculating the time, it is assumed that the magnetic-tape-search (mts) opera¬ 
tion (but not tape read and write) may proceed in parallel with other computer operations 
until either the search is completed or another tape instruction is given. The time for mul¬ 
tiplication and division, when the address section of these instructions refers to the perman¬ 
ent table of constants, ^ is made proportional to the power of ten involved and in all cases is 
less than the time required to perform the same operation using an ordinary SAC location. 
This is done because such multiplication or division would ordinarily be accomplished us¬ 
ing a numerical-shift operation, which in many computers requires less time than straight¬ 
forward multiplication and division. In all other cases, the time is calculated assuming 
that no operations may proceed simultaneously and that the time required is independent of 
the numerical quantities and memory registers employed. 


1. Actually, this integer is: (contents of excess) x 2^ 7 ♦ (contents of accumulator). 

2. See above. 
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Automatic Post Mortcms 

The SAC conversion and computation post mortems are virtually identical to those 
of the SS computer. The computation post mortem in SAC indicates the positions of the four 
tape units and specifies the most recently selected unit, if the tape has been used. No other 
changes in the post mortem have been made necessary by the modification of the computer 
instructions. 

197 THREE-ADDRESS COMPUTER 

A hypothetical three-address computer (called TAC) was developed during the past 
quarter by staff members of the S&EC Group. TAC was used during part of the 1954 sum¬ 
mer-session course 6.531 (Digital Computers: Business Applications). A description of 
TAC is given in Section 5. I of this report. 

198 STUDENT PROBLEMS CODED FOR SAC AND TAC 

This problem number was set up for the use of students enrolled in the special sum¬ 
mer-session course 6.531 (Digital Computers; Business Applications). 

All of the students in the course coded and debugged a problem assigned for solu¬ 
tion on the simulated three-address computer called TAC (developed under Problem 197). 
The problem consisted of determining and printing out the day of the week corresponding to 
an arbitrarily specified date. 

Most of the students also worked on a second problem to be solved on the simulated 
single-address computer, SAC (see Problem 196). For this second problem, the students 
were permitted to select from among three suggested problems covering the maintenance 
of bank-balance records, payroll calculations, and an autocorrelation calculation to detect 
a seasonal sale variation. 

199 LAMINAR BOUNDARY LAYER OF A STEADY. COMPRESSIBLE FLOW 

IN THE ENTRA NCE REGION OF A TU BE 

In connection with research on heat-transfer coefficients, recovery factors, and 
fraction coefficients for supersonic flow of air in a tube, a theoretical investigation of the 
characteristics of the laminar boundary layer in the entrance region has been carried out 
by Dr. T. Y. Toong of the MIT Mechanical Engineering Department. Partial differential 
equations of continuity, momentum, and energy were developed for the boundary layer. 
These were then transformed into a series of ordinary differential equations to be solved 
by WW1 for several entrance Mach numbers and thermal conditions at the tube wall. 


First, solutions are to be obtained for the case of constant viscosity and thermal 
conductivity. Then, the effects of temperature dependence of these prop*, sties arc to be 
studied. 

Solutions of the first set of differential equations for the case of constant viscosity 
and thermal conductivity have been obtained for six cases ol different entrance Mark num¬ 
bers and thermal conditions at the tube wall. This represents about 6% of the entire job to 
be done by WWl. 

Work is being done tonolve the second set of differential equations cor responding U> 
the same six cases. 

The coding has been carried out in the "algebraic system" developed by Dr. J. H. 
Laning under SfcEC Problem 108. This system makes use of Gill's modified fourth-order 
Runge-Kutta method. 

200 A STUDY OF RECURRENT EVENTS 

A study of recurrent events in sequences of binary numbers is being made by Miss 
Bergittc Jensen and Dr. G. P. Dinneen oi Lincoln Laboratory. In each experiment two 
Bcrnouilli sequences arc generated using a pseudorandom number generator. The basic 
sequence has a probability, p = P^. of a one which is small, less than 0. 1. This is denoted 
as a noise sequence. The signal sequence has a probability, p = P s » of a one which varies 
from about 0.2 to 0.9. The signal sequence is generated for a time which is short com¬ 
pared to the length of the noise sequence. 

Various sampling techniques are being tried to determine the optimum statistical 
test for determining the presence of the signal sequence and its position. The more com¬ 
mon statistical problem is that of deciding whether a finite sample comes from a noise re¬ 
gion or a signal region. Here we observe a finite sample, part of which is from noise and 
part from signal, and must determine the point of transition. Two methods are being tried. 
One is the success-run technique, where a particular sequence such as 0 0 0 0 0 is searched 
for. The other is the density method, where the density over some interval is measured as 
the interval moves through the sequence. Whenever the density exceeds a certain thresh¬ 
old, we say the sample is from a signal region. 

Programs for both of these methods have been coded and preliminary results ob¬ 
tained. It is expected that this study may be. completed during the next quarter. 

201 STUDY OF THE AMMONIA MOLECULE 

A self-consistent calculation of the electronic wave functions and binding energy of 
the ammonia molecule has been carried out. The method, which has been described in de¬ 
tail in previous Summary Reports under Problem 144, is this, each 1-elytron molecular 
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orbital is expressed as a line** ituu ui ihe »ei 01 atomic orbital* belonging to the 

nitrogen atoms and to the three hydrogen atoms. The total molecular wave function is the 
determinant with the 1-electron molecular orbital* a* elements. The expectation value of 
the Hamiltonian is expressible as a quartic form in the coefficients of the linear combina¬ 
tions. Insistence that the energy be stationary with respect to variation of these coefficients 
demands the diagonalization of an effective Hamiltonian matrix which it self contains the un¬ 
known coefficients. The problem is solved iteratively m what is termed a self-consistent 
manner. A first guess is taken for the coefficients, a new set is determined by the dlago- 
nalization of the effective Hamiltonian, this new set is fed back in, and so on until the input 
and output of a cycle agree to sufficient accuracy. 

The matrix components of the effective Hamiltonian a re composed of kinetic and po¬ 
tential energy integrals over the fixed set of 1-electron functions. Actually, these were 
taken to be not the primitive atomic set but or thogonalized and symmetrized combinations 
of them. The integrals over the atomic functions were evaluated on desk computers and 
IBM equipment. The integrals over the symmetrized combinations were obtained by means 
of a congruent transformation applied loa matrix whose elements were the atomic integrals. 
This was done on Whirlwind (a 16 x 36 congruent transformation). 

With these integrals, the self-consistent calculation was done on Whirlwind. Var¬ 
ious criterions of self-consistency were tried so as to judge the stationary property of the 
energy. The results are discussed by Dr. H. Kaplan in the 15 September Progress Report 
of the Solid Slate and Molecular Theory Group at MIT. 

202 CALCULATION OF VERTICAL ANTENNA COVERAGE SKELETON 

For several months investigations on vertical coverage interference patterns have 
been carried on by A. F. Barlholomay of Lincoln Laboratory. Complete detailed discussion 
on this subject can be found in his forthcoming report entitled "On the Calculations of the 
Vertical Coverage Interference Pattern." 

In brief, however, an interference pattern can be analyzed into two components: (1) 
the skeleton structure which consists of axes through the center of all lobes and loci of all 
null points forconstant point differences, and (2) eithera family of envelopes of cor respond¬ 
ing lobe points (the general case); or (in the flat-earth case) sinusoidal segments approxi¬ 
mating the actual lobe shapes. The difficulty in obtaining one complete pattern, notwith¬ 
standing several patterns which actually are needed, is twofold.(1) numerous formulas are 
required to determine a single point; and (2) the number of lobes in such a pattern is di¬ 
rectly proportional to the height of the antenna (hj) above the surface. For example, for an 
antenna height of 10,000 ft and a 100-mc radar, there will be roughly 2000 lobes. If we al¬ 
low 12 points as sufficient to describe each lobe, we would nerH £1,000 sepaiale calcula¬ 


tions for the complete diagram. 

Because of the lengthof the computations, programming the problem for Whirlwind 1 
seemed the best approach. Since the entire problem i* time consuming even for a large- 
scale digital computer, we considered only obtaining the skeleton structure of the pattern. 
We Hid develop a successful program in Uie (24,6) system where, tor a given antenna height 
of 10,000 ft, a skeleton pattern of the first 30 lobe lines was obtained. More explicitly. the 
program printed out in decimal form the and Y q coordinates for each n lobe line (1 ^ n 
^ SO) and for the curve of the earth's surface; and, at the same time, using the automatic 
curve-plotting routines, the program gave a scope display of the entire skelrton structure 
on one frame. The actual computer time consumed less than 10 minutes; equivalent cal¬ 
culations performed manually took several weeks. However, unforeseen complications in 
the original method of computation made us question the continuation of the problem on 
Whirlwind I. For the program to operate, limits must be set on dj, which is a iunction of 
n but which varies greatly for different n. Only through performing manually the entire 
calculation can such limits be found. Therefore, because of this and other inherent diffi¬ 
culties in the problem itself (sec above report), we conclude that, for the present, this 
problem is best handled by desk calculations. 

204 EXCHANGE INTERVALS BETWEEN REAL SLATER ORBITALS 

It is reasonable toexpect that our ideasof chemical behavior can be sharpened when 
we are able to compute this behavior from fundamental theory. A theory that ascribes cou¬ 
lomb forces, spins, and correct statistical properties to atomic nuclei and electrons with 
quantum mechanics in the form of Schroedinger's equation is believed to be sufficient to ex¬ 
plain all important chemical phenomena. Cases have been computed (Hylleras - the helium 
atom; James and Cooleridgc - the hydrogen molecule) to demonstrate the remarkable agree- 
ment of deduced behavior with expcrimcntal phenomena. Nevertheless, todate no ambitious 
computations closely related to the fundamental theoretical basis have been essayed. 

With simplifying assumptions, it is felt that such computations can be made feasible 
if the distribution (or wave functions) of the electrons are expanded in quickly convergent 
set of functions which form a complete set and, as orbitals of a single atom, have in a sin¬ 
gle function almost uniquely the properties of an orbital electron. Slater orbitals are such 
functions. 

An initial and computationally most difficult step in such a procedure is to compute 
integrals where the integrand is the product of two or four such orbitals and an operator. 
Il is this step that gives impetus to the "Integrals Project" now being carried on at the Uni¬ 
versity of Chicago under professors Mulliken and Roothaan. It is the eventual aim of the 
project to have all integrals necessary to molecular computations readily accessible. 
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The specific integral to be computed on Whirlwind I is a two-center exchange inte¬ 
gral, which is most easily explained in terms of the equivalent electrostatic problem. 

This problem involves the computation of the electrostatic interaction of two charge 
densities, fl(l), fl (2). by Coulomb's law. 


a (QflU) 

r i2 


dV ,dV . 
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m ■$ 0 


where the volume integrations are over all space. Each charge distribution is the product 
of two functions of the form 
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where the parameter ranges arc / f 1 ^ n. n = /, and ( ^ 0. The two functions whose 
product makes a charge distribution are located about distinctly separate origins, 

A transformation of the variables to elliptic coordinates is effected and— L_ ex¬ 
panded by the Neumann expansion. r li 

The resulting integral can be separated allowing the construction of an infinite series 
of integrals where the quadrature has been reduced from a sixiold to a double quadrature, 
the limit of the first quadrature being the argument of the second. 
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a is a function of thef's arising fromil(l) and R. a is similarly defined for £1(2). 
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where w^B.x) is a polynomial in integral powers of x. the coefficients of which arc de¬ 
termined as functions of all the parameters arising in a given charge distribution. 

The coefficients arising in the w /(B.x) polynomials are constructed from a set of 
spherical Bessel functions of imaginary argument and 1/2-int'ger order. Two arguments 


will be necessary for each total integration, and the process is to be carried out only once 
during the computation of a total integral. These functions are derived from the initial con¬ 
struction of four (unctions by means of a recursion formula. From these Bessel functions 
further functions will be constructed with the construction of each Iy by means of a simple 
recursion. At the outset it will be necessary to construct a table of coefficients. This ta¬ 
ble is repre sentative of the multiplication of the two \p functions together. The major portion 
of the computation is the numerical quadrature necessary to obtain .in 1^. The integration 
from 1 to cd is converted to an integration from 1 to 0 to give a finite range to the integra¬ 
tion. The procedure will be to evaluate the integrand (including the inner integrations) ex¬ 
actly and approximate the outer integration by some "numerical quadrature" procedure. 

To date, the initial coefficient table -- to be used in constructing the Wy^(B t x) pol¬ 
ynomial -- has been coded by P. Merryman and is now being tested. 

In view of the fact that the heaviest portion of the computational labor is in the nu¬ 
merical integration, rapidly converging numerical-quadrature schemes have been investi¬ 
gated. Of these schemes the most powerful, Gaussian, has proved adaptable, mutalus 
mutandis, and gives promise of decreasing the total necessary computations by the order 
of 70% of those necessary when using Simpson's rule. This scheme has been flow dia¬ 
grammed for the specific case and will shortly be coded and tested. 

A subroutinr giving spherical Bessel functions of imaginary argument and l/2-inte- 
gc r order, as well as square roots and exponentials, in (24,6) has been completed and tested 
by F. J. Corbato of the MIT Physics Department and is now ready for incorporation into 
this scheme. 

The procedure for computing an exchange integral is now sketched in sufficient de¬ 
tail. The remaining effort will be almost entirely devoted to the tasks of coding, testing, 
and assembly. 

Reference: Ruedenberg, Klaus, "The Two-Center Exchange Integrals" (Part II of "A Study 
of Two-Center Integrals Useful in Calculations on Molecular Structure"), Jour¬ 
nal of Chemical Physics, Vol. 19. No. 12, pp 1459-1477, December 1951. 

205 ELECTRON LATTICE INTERACTION IN SOLIDS 

Because of the large dipole moments associated with the long wavelength optical 
modes of vibration in ionic crystals, the coupling between these modes of vibration and the 
motion of a conduction-band electron is strong, too strong to be treated by perturbation the- 
ory. It was even conjectured at one time that such an electron could be trapped in the po¬ 
tential well created by the polarization charge it itself had induced. Although no experimental 
evidence for such a self-trapped electron has been found, it is still expected tha* the effective 
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mass of a conduction electron together with its associated lattice polarization could be very 
large. 

Various approaches to the problem of a single electron interacting strongly with the 
optical modes have been proposed. Frohlich, Pelzer, and Zienau* and Cross'^ have worked 
in a Fock space for the phonon assembly introducing cutoffs above 1. 2. or 3 quanta, then 
attempting an exact solution of the approximate problem remaining. Lee, Low, and Pines * 
and Gurari* have applied intermediate coupling theory in what amounts to a Hartree-Fock 
approach in the momentum representation of Fock space. Various Russian writers, includ¬ 
ing Bogolyubov, ^ Tyablikov. h and Pekar, ^ have worked in an electron-lattice configuration 
space making various approaches to an adiabatic approximation. 

Basic to all these approaches are certain assumptions. The periodic potential of the 
lattice is taken into account through the introduction of an effective mass. The lattice is 
treated as a harmonically vibrating continuum of pola rization charge in a phenomenological 

„ g 

Hamiltonian introduced by Frohlich. The analyses are restricted to small momenta, i.e., 
slow electrons. 

Q 

Several questions have been raised. Haken in treating a highly specialized model 
has suggested that the effect of the periodic potential may be considerable. Lee and Pines 
claim to obtain an exact solution in the limit of strong coupling, yet Pckar, using a varia¬ 
tional method, seems to obtain lower energies for sufficiently strong coupling. Lastly, the 
Fock-space methods suggest an effective mass increasing sharply with momentum, but these 
method* are valid only for small momenta. 

hi order to resolve some of these difficulties it is desired to construct a model which 
can be solved with sufficient exactness to serve as a yardstick in determining the validity 
of the various approximations. A possible model has been investigated exactly by Gross, ^ 
who determined the solutions to the eigenvalue problem arising when an otherwise free elec¬ 
tron is allowed to interact with only one traveling lattice wave and all motion is confined to 
one dimension. A possible limitation of this model is its spatial asymmetry. For suffi¬ 
ciently strong coupling, for example. Gross finds the state with minimum energy has non¬ 
zero momentum. Conceivably this is a fault of the model's asymmetry. To determine 
whether or not the asymmetry is truly a limitation to Gross's model, the problem of a free 
electron interacting with two degenerate modes (right- and left-going) in one dimension has 
been set up by T. D. Schulz of the MIT Solid State and Molecular Theory Group. The ma¬ 
trices of order arising on neglect of all states having more than n phonons ex¬ 

cited are being diagonalized for divers values of coupling, oscillator wavelengths, and total 


* Numbered superscripts refer to references at the end of this problem. 
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momentum (using a routine developed by Dr. Mvckler under Problem 134) with the follow¬ 
ing aims: 

1. To see how the effect of the inte raction on the ground state depends on wavelength; 
Z. To see if the asymmetry of Gross's model is important; 

3. To provide an accurate yardstick for evaluating adiabatic methods used on real 
crystals in the event the Gross's model is invalidated; and 

4. To provide a starting point for considering the effects of a periodic potential and 
a discreteness in the vibrating lattice. 

Two values of the coupling constant (intermediate and strong) and two wavelength* 
(intermediate and long) are being investigated. For intermediate coupling good convergence 
has been obtained. For strong coupling it is found that a better criterion than total* number 
of quanta must be used in cutting off the infinite secular equation to obtain good convergence. 
Dr. Meckler is considering modifications in the present routine which would include only 
those states with n lowest diagonal matrix components. 
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206 ELECTRON ENERGIES OF THE MOLECULE H^ 

The complete problem, undertaken by A. Dalgarnoof Queen's University of Northern 
Ireland in association with the MIT Physics Department, was the evaluation of the rates of 
the collision processes 
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(1) charge transfer 


(2) associative detachment 


/ -flr. ^r,. 

/* ** 0/ r u > dr, tir 2 


and, in addition, (3) the evaluation of the mobility of H tn .» gas nf atomic hydrogen. Such 
processes are of importance in the physics of the upper atmosphere and of the sun. 

The basic data required for (1), (2), and (3) are the electronic energies of the two 
lowest states of the molecule H , . (It may be mentioned incidentally that the calculations 
will also provide useful information on the dissociation energy and electron affinity of H^~.) 
The energies of H^ can be expressed in terms of certain molecular integrals, and it is for 
the evaluation of these that the Whirlwind I computer has been used. 

The integrals involved have the forms: 

/ n /"a , n -0 r b . f n ‘"'a l n ^ T b A 

/ r a r b e dr , / r e V r.e dr . 


f n - ar la n'0 r lb k 

y r la e r lb * r Za e r i 


’(l/r u ) dr, dr 2 


< >/ r 12 ) dr, dr^ x 


By an elementary expansion of 1, in terms of spherical harmonics, it may easily be 

shown that all these integrals are expressible as somewhat lengthy combinations of prod¬ 
ucts of auxiliary integrals 


A ( X) = 


«,<,>■ f 


where r, (or r ), etc., is the distance of electron I from nucleus A, r,, is the distance 
ia a 12 

of electron 1 from electron 2, a . £ , y . & arc parameters, and m, n, k, i are integers 
which may be negative. Depending on the values of the various parameters, essentially six 
different integrals arise, known as 

1. Overlap 

2. Coulomb attraction 

3. Kinetic energy 

•4. Coulomb repulsion 

5. Hybrid 

6. Exchange 

and they are required for a Urge numbei of values of a , £, y, B and the nuclear separa¬ 
tion A of which the integrals are a function. 

Integrals 1, 2, and 3 were already available in a program written by F. J. Corbato 
of the MIT Physics Department. It was hoped to program 4, 5, and 6, but because of time 
limitations only 4 and 5 were completed. However, 6 is being programmed under Problem 
204 by P. Merryman. 

The integrals which were programmed and for which results were obtained are spe¬ 
cifically 


A routine for generating these functions was available. 

208 INTERCEPTOR FLIGHT-CONTROL PROBLEM 

A system of 14 first-order nonlinear differential equations representing an inter¬ 
ceptor flight-control problem was solved by M. Merwin of the MIT Dynamic Analysis and 
Control Laboratory using the Gill adaptation of the fourth-order Runge-Kutta method. The 
Whirlwind solution showed that the equations were unstable. These results, reproduced 
later on an analog computer, proved to be most useful, since the system defined by the dif¬ 
ferential equations was thought to be an optimum system. 

209 NUMERICAL SOLUTION OF HOMOGENEOUS LINEAR DIFFERENTIAL 
EQUATIONS IN THE QUADRATIC POLYNOMIAL COEFFICIENTS 

This problem was selected by Dr. J. C. P. Miller of Cambridge University, Cam¬ 
bridge, England, as one that not only would provide a generally useful and interesting rou¬ 
tine but could be completed during his brief visit at the Digital Computer Laboratory. 
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The equation to be solved is 

♦ a>) i-* ♦ (b l » Z tb 2 x +bj)^ ♦{c 1 , i t c J «tc } )y.O, 
d v 

starting at x = x Q , with y = y Q . = y^*, at interval h in x. 

The method of solution is to obtain a recurrence relation forterms in the Taylor ex¬ 
pansion at each x and to find values of y and y' at x + h, by summation, and then start again 
with the new value of x. 

Writing the equation as 

p(*> y z * q(*) y, + r(*) y = o 

and differentiating n times we get 

py n ♦ 2 + < np ' * q > ♦ 1 f P " + " q ' f r) y n + q” ♦ "r')y n . , 


+ =inj_L) r .. y 


y 


Then if we write T = -" or y = n'. T h' n , we obtain 

n nj 'n n 

(n r l)(n + 2)pT R + 2 + (n + l)(np' *t)hT n + , + O p.. + nq . f r)h Z T n 

♦ (— 5 — q" + r')h 3 T . + 1/2 r"h*T , = 0. 
c n - i n * t 

This defines T n + 2 for n * 0. I. 2. ... if we take T 2 - T j = 0, T Q = y. and Tj = hj. 
Then 


y(x th) = T q (x + h) * T„ + Tj + T 2 + ... 


hy' (x + h) = T|(x + h) = T, + Z1 z + 3T } + ... 


gives the new T fl and Tj. The processes cease when two successive terms are negligible. 
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This equation covers many familiar and unfamiliar functions, it has been used with 

2 

e satisfying y - 2xy‘ - 2y = 0 

cos x satisfying y" ♦ y - 0 

/* 

J c dx satisfying y" - 2xy' = 0 
J o 

x 1(J satisfying xy" - 9y' » 0 

y = log x satisfying xy" + y' - 0 

showing that nearly seven figures can be readily obtained. The greatest error occurred 
with e X , where units in the seventh figure accumulated in eight steps, corresponding to a 
factor exceeding 10^ in the increase of the function. 

211 SERVO RESPONSE TO A COSINE PULSE 

Th*« problem is concerned with extending the curves in Figure 25-16, page 705 of 
Instrument Engineering, Volume 11. by Draper, McKay, and Lees. It is desired to find the 
maximum response x(t) for 1440 parameter pairs (£, T^/ T n ) in the equation 

(1/w “) ‘j l -f (2f/w ) x f x = G(l), where T = 2f/w , 

X = dx/dt. * = d 2 x, dt 3 , G(t) = (1/2) [j - cos (2irt/T p )] for 0 5 I 5 T , 

and G(t) = 0 for t > T^, 

The parameters £ and T T are from the range 0.01 ^ £ < 5, 0. 1 £ T /T ^ 10. 
r 3 p n 9 p n 

Analytic expressions for x(t) may be found in the above-mentioned book, but these 
were found *mprat.lical for determining the maximum response. Accordingly, a step-by- 
step approximation to the response curve using Taylor scries is employed. The step-by- 
step procedure is stopped when the derivative becomes negative, and an interpolation is 
carried out back into the last step. This procedure was applied using parameters different 
from those of the present problem by Dr. J. H. Laning, Jr., in coding the problem for a 
card-programmed calculator. Coding of the problem using the CS II system was accom¬ 
plished by Dr. J. M. Stark of the MIT Instrumentation Laboratory. 
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3.1 Systems Engineering 


During the past quarterly period the Whirlwind Computer System reliability was 
high. According tv calimales by computer operators, 96% ot the hours ot assigned 

operation time was usable. A detailed account of the assignment of computer time to the 
various major activities is given in Fig. 3-1. 

Work by the systems group on the central computer has included modernization of 
the arithmetic-element control, reorganization of the power-supply control fora large seg¬ 
ment of the terminal equipment, and some equipment relocation and circuit changes for 
easier, less-frequent trouble shooting. Marginal-checking facilities, extended to cover 
more of the terminal equipment, have already proved their value by revealing some circuit 
shortcomings now corrected. 

Two additional drum groups of auxiliary storage added to the 15 already available 
to the computer increase drum storage from about 31,000 to about 3b, 000 registers. 



Hours Per Week 

Activity 

July ] August | September 

Total 

Hour* 


4-8 

15 

16-22 

21-49 

10-5 

6-12 

11-19 

40-26 

27-2 

1-9 

10-16 

17-2! 

44-30 

Marginal Checking 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

52 

Installation 

0 

4 

8 

8 

10 

8 

8 

9 

8 

0 

7 

7 

8 

86 

Maintenance 

15 

11 

22 

25 

30 

36 

27 

20 

22 

11 

38 

36 

26 

3 39 

Terminal Equipment 
Testing 

48 

43 

39 

3S 

28 

48 

30 

10 

14 

30 

29 

25 

33 

412 

technician Instruction 

5 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

53 

Scient.dc and Engineer¬ 
ing Computation 

4b 

47 

50 

50 

49 

41 

_ 

41 

60 

50 

43 

33 

— 

43 

41 

*91 

Other Applications 

25 

14 

21 

21 

24 


26 

22 

27 

19 

19 

2) 

31 

293 

Total Hours 

123 

148 

148 

147 

149 

142 

_ 

152 

149 

149 

96 

134 

142 

149 

1848 


Most of the circuitry and logic of the arithmetic element of the computer iscontem- 
porary with the original system installation. The elementary state of the art at that time 
led the origina 1 designers to build into their circuits a large degree of flexibility against 
unforeseen circumstances. Experience in computer operation has shown some of the orig¬ 
inal equipment to be unnecessarily complicated; where practicable WWI is being improved 
in this respect. 

Revisions of the arithmetic-element control included standardization of circuitry 
and chz;.ges in logic. These modifications to the stopcounter and divide control eliminated 
excess circuitry, reduced the time required for the divide operation by one-third, and in¬ 
creased the operating margins of the system. 

A large part of the terminal equipment for WWI, one floor below the main section 
of the computer, was installed at a later date than the bulk of the WWI circuitry, and a sep¬ 
arate power-supply-control system was designed for it. It has long been evident that the 
power-supply-control system left much to be desired because of the inaccessibility of the 
circuitry and because the power is not turned off in proper sequence in all modes of shut¬ 
down. This results in writing between the slots on the magnetic drums and has materially 
reduced drum reliability. In addition to these shortcomings, it is now felt that the present 
independent control of the auxiliary drum and bufler drum is not necessary. 

A new power-supply-control system eliminates these objections mentioned and re¬ 
duces the number of independent controls for the two groups of equipment. 

The magnetic-tape equipment in use with WWI can be logically divided into two sep¬ 
arate systems, the auxiliary memory (which has the exclusive use of two tape units) and the 
delayed-printout system (which shares the other three with the computer). The magnetic- 
tape system, in its auxiliary-memory capacity, has been operating satisfactorily with the 
computer during the past quarterly period. The delayed-printout equipment has been giving 
some trouble, and several changes were made in the system to improve its reliability. Cir¬ 
cuit modifications improved the operating margins, and a relay in the control register al¬ 
lowed the Flexowriters more time to operate. Additional relays have been installed to re¬ 
duce unwanted signal coupling. The equipment is being moved to a new position where the 
panels will be more accessible for servicing. 
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Fig. 3-1. Allocation of Computer Time 
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In the past two years the inpul-output by stem for WWi has grown so rapidly that mar¬ 
ginal-checking facilities were sometimes not added as new equipment went into service. 
When this omission was made good, marginal-checking lines frequently were not optimized 
in their arrangement. During this quarter new marginal-checking circuits were installed, 
and old lines for the in-out system were reassigned to bring them up to the standards of the 
rest of the computer. 

Some developmental work was necessary prior to this installation because of the na¬ 
ture of several of the circuits associated with the terminal equipment. A method of mar¬ 
ginal-checking the crystal gates in the "and" and "or" circuits for the intervention regis¬ 
ters was developed so that it is now possible to pinpoint weak crystals in these registers. 

When marginal checking was extended to the gas-tube indicator-light register, hith¬ 
erto unsuspected weaknesses of the circuit were revealed by the low and rapidly deterior¬ 
ating margins. As a result of this discovery, the extinguishing circuit was redesigned to 
make the tube requirements less critic'1 and to give more reliable operation. 

Auxiliary-storage capacity of the buffer drum has been increased from about 8000 
to about 12,000 registers, each of 17 bits, by using Groups 2 and 3 for this purpose. 

The source of inte rmittent c rosstalk between the buffer- and auxiliary-storage sec¬ 
tions of the buffer drum has been located and removed. Occasional buffer-drum parity 
alarms over a period of several weeks are attributed to this fault. The change should im¬ 
prove the reliability of the buffer drum to the point where it approaches that of the auxili¬ 
ary drum. 


4.1 Vacuum Tubes 

•4. 1. 1 Vacuum-Tube Life 

During the third quarter of 1954 the WWI computer operated for 1680 hours. 
Vacuum-tube life has been calculated for five different tube types as outlined in 
Summary Report No. 36. A summary of this information together with previous, data is 
shown below. 

FAILURE RATE. PER CENT PER 1000 HOURS 


'l ube Type 

1952 

1953 

First 

Qua rter 
1954 

Second 

Quarter 

1954 

Third 

Quarter 

1954 

7AD7/SR 1407/6 145 

2.00 

J. 3 

1.75 

1.6 

1.4 

7AK7 

0. 26 

0. 43 

0. 5 

0. 6 

0. 3 

6080/6080WA/. A 


6. 6* 

1. 1 

1.5 

1.4 

5965 


0.2» 

0.4 

0. 3 

0. 2 

6BL7GT 


0.7* 

0. 3 

0. 5 

0. 5 


♦Last quarter 1993 only 

There are no major changes in these figures from the second quarter. 

Several life tests are still active. That of the GE 5969 has reached 9400 hours with 
one section conducting and the other cut off. Monotonic shifts in contact potential, which 
have caused the grids of the cutoff sections to become more positive and those of the con¬ 
ducting sections to become more negative, have occurred. The average difference between 
the sections is about 0.7 volt. Considerable grio emission has developed on the cutoff sec¬ 
tions while the grid currents of the conducting sections are normal. No inierface impedance 
has been observed. 

A life lesion the RCA A4688B (developmental 9969), being run with one section con¬ 
ducting and the other cut off. has reached 2000 nours. Some grid-cathode leakage has be¬ 
gun to develop, predominantly on the conducting sections. Four of the 19 cutoff sections 
show interface impedance ranging from 1 to 3 ohms with a fifth section showing 16 ohms. 
No interface is present on the conducting sections. 

The life test of special 5687's with Ail Cathaloy (tungsten-nickel) cathode sleeves 
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has reached 3400 hours with one section conducting and the other cut off. Considerable grid 
emission has developed on the cutoff sections; grid currents in the conducting sections are 
low. The difference in plate curren. at zero bias between the conducting andcutoff sections 
of a given tube is quite small. This is a considerable improvement over the standard 5687. 
Three of the ten conducting sections show interface impedance ranging from 3 to 15 ohms, 
while no interface is apparent on the cutoff sections. 

A second life test involving tungsten-nickel A3 1 Cathaloy cat »ode sleeves, being run 
on some special GE Z-2177's (improved 5965), has reached £000 hours with one section 
conducting and the other rut off. In general the grid currents are low indicating that grid 
emission may not be a problem. Interface impedance ranging to 55 ohms is present in £1 
of the 35 conduction sections, the cutoff sections show a maximum of 15 ohms in five of the 
35 sections. Worthy of note is the fact that there is more interface impedance present on 
the conducting sections of this tungsten-nickel alloy cathode than on tne cutoff sections, as 
is also the case on the 5687 tungsten-nickel alloy samples. This is just the reverse of ex- 
peri* nee on active alloys where silicon is the impurity. 

4. 1.2 Vacuum-Tube Research 

An investigation of 30,000-hour 7AK7's from the five-digit multiplier has been 
completed and a report, M-3020, published. Considerable interelectrode leakage probably 
caused by sublimation on the micas from the cathode nickel sleeve was found. On tapping, 
transient decreases in intcrelectrodc resistance probably caused by changing contact of the 
side rods with the micas resulted in more than half of the 30 tubes indicating intermittent 
shorts. The sensitivity of the detecting equipment is such that a maximum resistance of 
approximately 5 megohms lasting at least 500 microseconds can be detected Barring ca¬ 
thode problems, interelectrodc leakage may be a serious limitation on long tube life. 

An extensive investigation was made of the state of the cathodes of these tubes. To 
eliminate effects of geometry, the tubes were triode connected to +90 volts and the control 
grid pulsed at 1% duty factor up to +30 volts. A useful criterion of the state of the cathode 
which evolved from the investigation was the change in pulsed triode-connected plate cur¬ 
rent at a fixed high value of grid drive ( + 25 volts) while varying the heater voltage from 
6.3 to 5.0 volts. This value ia called "n" and is given as a percentage of the 6. 3-v cur¬ 
rent. From the basic assumption that an increase in n is caused by an increase in the over¬ 
all work function of the cathode, it is possible to plot a family of curves showing variation 
in the triodc-connccted plate current as n varies; new tubes show an n of 6%. 

By adding to the measured current of the old tubes the percentage decrease calcul¬ 
ated from their n measurement, a distribution of currents is found whose mean is about 5 fl fc> 
below that of new tubes and whose standard deviation is essentially the same as that of new 
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tubes. A Ficher "t" test, however, shows a fair probability that the 5% diffe rence could be 
attributed to chance. Now by comparing the percent decrease from a perveance curve at a 
fixed triode-connected plate current (saturation) with the values predicted from the n cal¬ 
culation on several representative tubes, a correlation coefficient of 0.97 is found, indica¬ 
ting a rather high degree of certainty for the validity ol the theoretical drop. 

The n test has proven itself to be a rather valuable tool in our vacuum-tube work in 
d-termining whether ? given low-current tube is a victim of a poisoned cathode or some 
other phenomenon. 


4.2 Component Replacements 

Fig. 4-3 lists the replacements of components other than lubes during the third 
quarter of 1954. 
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Component 

Type 

I oial in 
Service 

No. of 
Failures 

Hours of 
Operation 

Comments 

Capacitor 

Electrolytic 
4.0-pl 450-v 

136 

1 

4000 - 5000 

Low capacitance 

Crystal 

Diodes 

IN34A 
(glass body) 

16.881 

2 

2 

1 

2 

21 

2000 - 3000 

5000 - 6000 

6000 - 7000 

7000 - 6000 

8000 - 9000 

Low R. 

Low 

Low R“ 

Low R? 

Low R b 

D- 358 

(ceramic body) 

5137 

1 

25000 - 26000 

Low 

IN92 

diffused junction 

61 

1 

3000 - 4000 

Drift 

Rectifier 

Selenium, 

IRC. 19PA6. 
156-v d-c. S-ma 

SM 

2 

4000 - 5000 

e 

B 

a. 

O 

Relay 

Allied 

BODO-bD-38 

675 

l 

1000 - 2000 

Failed to 
operate 

Resistor 

Deposited carbon. 
9000-ohm 
1/2-watt fl% 

878 

1 

1 

0 - 1000 
2000 - 3000 

Above tolerance 
Above tolerance 

Switch 

Toggle. SPST 

606 

1 

24000 - 25000 

Intermittent 

contact 

Trans - 
former 

Pulse. 

Hipersil core, 3:1 

4280 

‘ 

2000 - 3000 

Open primary 


Fig. 4-2. WWI Component Failures July 1 - September 30, 1954 
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5.1 MIT Courses 

Summer Session, 1954 

During the summer of 1954, MIT offered 35 different special summer programs, 
courses of l or 2 weeks duration largely for the benefit of interested people from business 
and industry. Two of these were directly concerned with digital computation, both being 
arranged by Professors Charles W. Adams and Stanley Gill with the full cooperation of 
the Digital Computer Laboratory. Several of the other programs were also of considerable 
interest in the field of computer application. For example, there were courses in Analogue 
Compulation, Numerical Control of Machine Tools, Operations Research, and Mathematical 
Problems of Communications Theory. 

The first of the digital-computer subjects, a 1-week seminar entitled Digital Com¬ 
puters: Advanced Coding Techniques, was intended as a meeting ground for experienced 
computer programmers interested in summarizing the present state of the art and in ex¬ 
changing ideas for future work in automatic coding. Planned for 25 persons, the course 
was with great difficulty limited to 42. including many experts in the field. While concrete 
results are hard to summarize, it was felt to be a fruitful experience for all concerned. 

An introductory 2-week course in Digital Computers: Business Applications stressed 
coding techniques and other aspects of computers in the light of commercial and industrial 
management information handling. The 61 students divided into four broad clarses, with 
25 representing 16 varied businesses. 17 representing 7 insurance companies, 7 repre¬ 
senting a number of government agencies, and 12 representing 6 computer manufacturers 
and other special interests. To facilitate the teaching and use of practical computers 
typical of machines commercially available for business use, the summer-session com¬ 
puter developed for the 1953 summer session was revamped into a single-address com¬ 
puter (SAC), andan entirely new three-address computer (actually a 3 4 1 address machine), 
called TAC, was developed. 

These two hypothetical computers proved very successful in demonstrating the 
v»rinn« features of digital computers, with only a small attention to the many trivial, 
almost unavoidable,details that tend to clutter up any practical computer design. 

Fall Term, 1954 

The principal course on machine computation being offered at MIT in the fall of 
1954 is 6.25, Machine-Aided Analysis, a survey of computing techniques aimed largely at 
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seniors in Electrical Engineering. This subject, first offered in the spring of 1954, had a 
fall-term enrollment of about 55 seniors and graduate students. It consequently was divided 
into three 20 - man sections taught separately by Professors Lmvill, Booton, and Adams. 
Practice problems were planned to allow each student to use both a REAC and the Whirl¬ 
wind 1 computer (simulating the hypothetical I AC which had been developed for the 1954 
summer session). Exercises using desk calculators and a card-programmed calculator, 
which were included in the first presentation of the subject, were eliminated in this second 
presentation to permit more time for studying techniques of problem analysis and of error 
analysis. 


5.2 TAC--A Three-Address Computer 

TAC is a hypothetical computer developed for use in MIT's special summer pro¬ 
gram in Digital Computers: Business Applications, August 1954. A compiler-interpreter 
program for MIT's Whirlwind I permits Whirlwind to simulate TAC. 

TAC Definitions 

A character is any of the following letters, digits, superscripts, symbols, or 
machine functions comprising the complete vocabulary of the MIT Flcxowriters (augmented 
by $ and * ). 

TC (('/)-+. 2 0 *468 0i468 ab...z 
RBH$:(_=,| S 13579 13579 AB. ..Z 

where R = return carriage, T = tab (stops every tenth space), B = back space, C = color 
shift (starts black), H - stop, I = ignore, and S = space. 

A word consists of nine characters, which are identified from left to right by the 
letters abedefghi or by the digits 1 to 9. 

A number is a word in which the first character, a, is -, -f, or 0, and each of the 
remaining eight characters bedefghi is one of nine digits 0, l, 2, ... , 9, forming an 8-digit 
decimal integer. 

An instruction is a word in which the first character, a, represents one of the 15 
operations listed below and the remaining characters, paired, form four "addresses." 
The scheme is shown below: 
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Operation character 

Address of first operand (except 
in N, L. R. P, T) 4. I 

Address of second operand- 

(except H. I. T) 


rr 

T 1 

- 




Address of next instruction 
(except in C, K, 1) 

.Address to which result is 
sent (except in C, K, P, R. 1) 


An address normally refers to one of 110 registers which are numbered 00, 01, 

. . . , 99, xO, xl.x9. However, there are certain exceptions to this: 

1. The be address associated with operations N, L, R, P, and T is omitted, and 
the charade rs are usedas a continuation of the operation designation (c. g. . N + , N-, below), 

2. The fg address for operation P is omitted, and the characters are used to con¬ 
tinue operation specification: 

S. The registers xl and xZ can be treated as a single register for certain purposes, 
and this register is denoted by the address xx. For operations A, S, M. D, or N, the 
contents of xx will be treated as a single 16-digit number, provided that xl and xZ both 
contain 8-digit numbers having the same sign. Thus if xl contains 412345678 and xZ con¬ 
tains +98765432, then xx will be deemed to contain +1Z345t>7898765432. 

We note that address xO always contains exactly zero. A post mortem (sec below) 
will always occur if any instruction contains an address which is not legitimate according 
to the above rules. 

A post mortem is performed automatically whenever TAC encounters any impossible 
instruction (in particular, any wffrd not starting with one of 15 legitimate operation code 
letters, or containing an illegitimate address, or violating the special conditions listed on 
pages 66, 67. ara 68). The post mortem consists of printing the location of the illegitimate 
instruction itself, the contents of the registers it refers to, a sequence table listing the 
locations of the instructions performed just prior to the post mortem, and an altered-word 
table listing the contents of all of the registers whose contents have been altered during the 
program. 

The time required to perform a given instruction maybe calculated from the follow¬ 
ing facts: TAC storage consists of a magnetic drum with four groups of 25 words (0-Z4, 
25-49, 50-74, 75-99) revolving at 40 revolutions per second (1 millisecond per word); 
access to the x registers is one word time (no waiting time required): and TAC always 
follows the control sequence: it acquires the instruction, acquires the operand (s*), per- 


* If two references to the drum are required to acquire operands or store results, TAC 
searches for both simultaneously, so that the time required is the longer of the two access 
times computed independently. 
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forms the required operation (in the number of word times indicated in the following list), 
and stores the result (*•) (omitting any storage accesses where none is required, as with 
the be address of N. L. R, P, or T). 

Symbols Used 

abedefghi represents the nine characters of a TAC word. 

C( ) represents "contents of." Thus C(bc) represents "the contents of the register whose 
address is be" (where b and c, the second and third characters of the instruction, 
must each be either one of the decimal digits or the letter x). 

—represents "becomes the new contents of." Consequently C(bc) t C(de)-^fg should 
be read as "the contents of be plus the ontents of dr becomes the new contents of 
register fg, replacing whatever was in fg but not changing what is in be or de." 
xx represents TAC's double-length register made by pairing xl and xl. 

> represents "is greater than." < represents "is less than." 

10 C (read as "ten to the c") represents a one followed by c zeros (for example, 10** = 

100 , 000 ). 


TAC Instruction Code 


Name 

Code 

Letter 

Word Times 
/ Ope ration 

Function 

Post Mortem Will Occur 

Read 

R 

6/char. 

Read enough characters from 
punched tape to fill the positions 
numbered b thru c in register de, 
without changing the other digits 
of C(de). 

(ignores deletions, deals properly 
but reads all other characters incl 
lines, tabs, and carriage returns 
each read as three characters) 

If tape contains illegal 
characters, or if b=0, or 
if c = 0, or if b>c, or if 
de = xx or is not a legiti¬ 
mate address. 

with upper and lower case, 
uding back spaces, under- 
explicitly, and $ and i are 

Print 

P 

10/char. 

Print the characters in the posi¬ 
tions numbered b thru c in regis¬ 
ter de, preceded by the character 
f and followed by g. 

Unless b = 1, 2, .... 9 and 
c = 1, 2. . . . , 9 and 
de is a legitimate address 
(not xx) or it b> c. 

Tape 

Read 

Tape 

Write 

Tr 

30/block 

Read ten words from block de (or 
from next consecutive block if 
de = 00) on tape unit c into regis¬ 
ters fg, fg f 1.fg + 9. 

Unless b = w or r 

" c = 1, 2. 3, 4 and 

" fg = 00, 01, ..., 

90 or xO 

" dc = 00, 01. 

or 40 

if de = 00 and block 40 has 
been just read or written. 

Tw 


Write ten words onto block de (or 
onto next consecutive block if de 
- 00) on unit c from registers fg, 
fg ♦ 1.fg 4 9. 
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Input 

I 

5/ char. 

Without altering the present con¬ 
tents of storage, start reading a 
punched tape containing a TAC 
program in conventional form, 
ending by taking the next instruc¬ 
tion from the address preceding 
"start" at the end of the tape. 


Halt 

H 


Stop computing. Start at hi only 
pressed. 

if the start button is de - 

Add 

A 

2 

C(bc) i C(de)-»fg 

1 . If the magnitude of 
the result exceeds 99,999, 
999 (or 9, 999, 999. 999, 
999,999 if fg xx). 

2. If be - xx or *de xx 
orfg xx and(be)and(de) 
are not both numbers. 

3. If the column-by-col¬ 
umn addition or subtrac¬ 
tion involves any nondigit 
N in any arithmetic oper¬ 
ation other than N f 0 - N, 
N - 0 = N, or N - N 0 
(see next page fo: details). 

Subtract 

S 

2 

C|bc) - C(dc)-»fg 

Multiply 

M 

10 

C(bc) x C(de)-»fg 

1. If magnitude of result 

19, 999 (or 
9.999,999.999,999.999 if 
fg = xx); 

2a. In M and D only, un¬ 
less C(bc) and C(de) are 
both numbers; 

Divide 

D 

25 

C(bc) ♦ C(de)-^fg 
(quotient rounded) 

Numerical 
Shift Left 

Shift Right 

N + 

? 

C(de; x I0 C —*.]g 

N- 

2 

C(de) ► 10 C —»fg 

2b. In N only, unless b = 

+ or -, c = decimal digit, 
and C(de) ^ a number. 

Compare 

Numerical 

C 

iy 

2 

Take next instruction from: fg if 
C(bc)>C(de). hi if C(bc)<C(de), 
or next register consecutively if 
C(bc) = C(de). 

If C(bc) and C(de) are 
not both numbers. 

Logical 

Shift 

L 

2 

Shift C(de) and C(fg) cyclically c 
places (left if b is +, right if b is 
-)• 

Unless b = f or c = 0, 

1.or 9; and de and fg 

arc legitimate addresses; 
r»r if he, de. or fg = xx. 


A post mortem occurs on A, S, M, D, C, K, and E unless be, de, and fg are all legiti¬ 
mate addresses. 
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Com pare 
Logically 

K 

2 

Take next instruction from fg if 
C(bc )> C(de), hi if C(bc) <C(de), 
or next register consecutively if 
C(bc) is identical with C(de). 

Here the symbol >is defined as fol 
actcrs of C(bc) with those of C(de 
the left end until identify is establis 
is found to be to the right of the ot 
which case the word containing sa 
greater than the other logically. 

If be, de, or fg = xx. 

lows. Compare the rhar- 
) column by column from 
hed or until one character 
her in the list below**, in 
id character is said to be 

Extract 

E 

2 

In those columns, and only those, 
in which C(de) has odd characters 
(listed in the lower line of the two 
lines below**), replace the char¬ 
acters of C(fg) by the characters 
occupying corresponding columns 
in C(bc), without altering the 
other characters in C(fg). 

* If be, de, or fg = xx. 


r~ 


_ . , „ 0 2 4 6 8 

n 


small 

even. 

T C t 


a b. . . z. 

large 

end 

odd. 

R B H 

$:( = .! S1S57<) 13579 AB..Z 

end 



— — 

_ - 

-J 



1. Register xO contains 000000000. If an instruction attempts to put other information 
into it, the information is lost. No post mortem occurs. 

2. The next instruction is taken from hi (unless otherwise specified in operations C and 
K) except in I. 

3. A post mortem always occurs if hi is not a legitimate address (except in I, and in H if 
not restarted). 

Addition and Subtraction in TAC 

The rules which TAC obeys in adding or subtracting are the following: 

1. While a positive number may start with zero instead of plus, a positive sum or 
difference formed by TAC always starts with a plus sign, 

2. If any or all of the addresses be, de, or fg in the instruction designate register 
xx, both C(bc) and C(de) must be numbers. If fg does not equal xx, the magnitude of the 
sum or difference must be less than 10®. 

3. If none of the addresses be, de, or fg designates register xx, the situation de¬ 
pends on whether addition or subtraction is required and on what the first characters of 
C(bc) and C(de) are. Summarized below are all possible combinations of such first char¬ 
acters The symbols used are: 

• represents either a plus sign or a zero 
4 represents a plus sign or a minus sign as appropriate 

k represents any character other than +, or 0 
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k represents any other character other than ♦, or 0 (different from k) 

... represents any set of eight characters to the right of the sign, i.e., in posi¬ 
tions bc2-9 or de2-9 

(w) - letters in parentheses refer to the first of three different cases tabulated at 
the bottom of the page, wherein the right-hand eight digits of the results are discussed in 
terms of C(bc2-9) and C(de2-9) 

P. M. indicates that no result is possible and that a post mortem always results in 
this situation 



if C(bc) = 

e... 

9. . . 

9. .. 

«... 

-... 


k. . . 

k. .. 

k. . . 

k. . . 

for 

and C(dc) = 

e... 


k. . . 

e... 

-... 

k... 

©. . . 

-... 

k... 

k'. . . 

ADD 

then C(fg) = 

»(w) 

+<*) 

k(w) 

♦to 

-(w) 

My) 

Mw) 

My) 

P. M. 

P. M. 


if C(bc) = 

o... 

e... 

e... 


-... 

. • 

k. . . 

k. .. 

k. .. 

k. . . 

for 

and C(de) = 

e... 


k. .. 

©... 

-... 

k. . . 

e... 

-.. . 

k. . . 

M.. 

SUBTRACT 

then C(fg) = 

+<*) 

+(») 

P.M. 

-<«) 

♦to 

P.M. 

My) 

M«) 

♦to 

P.M. 


The right-hand eight columns of the result, i. e., 


C(fgZ- 

9), are derived from C(bc2-9) and C(dc2-9) thus: 

A post mortem will occur 

W 

C(bc2-9)and C(de2-9) are simply added column 
by column. 

IF the sum exceeds 99,999.999 
or IF the addition requires add¬ 
ing a carry or any cha/acter 
other than zero to any nondigit. 

X 

If C(bc2-9) is larger than-C(dc2-9), assuming 
all nondigits to be larger than zero, then 
C(de2-9) is subtracted column by column from 
C(bc2-9), otherwise C(bc2-9) is subtracted 
column by column from C(de2-9). 

IF the subtraction requires sub¬ 
tracting a borrow or any char¬ 
acter other than zero or itself 
from any nondigit. 

V 

The quantity C(bc2-9), or C(dc2-9), whichever 
is preceded by the minus sign, is subtracted 
column by column from C(de2-9) or C(bc2-9). 

IF the difference so formed is 
negative or IF the condition of 
case x arises. 
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6. APPENDIX 


6.1 Publications 

Project Whirlwind technical reports and memorandums are routinely distributed to 
only a restricted group known to have a particular interest in the Project, and to ASTIA 
(Armed Services Technical Information Agency) Document Service Center. Knott Building. 
Dayton. Ohio. Requests for copies of individual reports should be made to ASTIA. 

The following is a list of memoranda published by the Scientific and Engineering 
Computations Group during the past quarter: 


No. 

Title 

Date 

Author 

SR-38 

Summary Report No. 38, Second 
Quarter 1954 



M- 289 1 

Proposed Change in WWI to Make 
Possible Programmed Recording 
of the "Lights" 

7-13-54 

A. Siegel 

DCL-7 

The "FL" Flexownter Code, Binary 
Numerical Sequence - reprint p. 20, 

M- 1623 



DCL-8 

The ”FL" Flexowriter Code, Alpha- 
numerical Sequence - reprint p. 19, 

M-1623 



DCL -1 1 

Automatic Scope Output Display 
Routines 

8-30-54 

N. J. Saber 

DCL-12 

A Program for Transferring Binary 
Information Back and Forth between 
Paper Tape and Magnetic Tape 

9-1-54 

S. F. Best 










